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Thesis entitled *Dilute and concentrated Mixed Valence 
Systems — Sane Simple Models r submitted by Keys Sur in 
partial fulfilment of the requirements of the Ph.D. 
degree to the Department of Physics, Indian Institute 
of Technology,- Kanpur. 

September 1978 

Some simple microscopic model Hamiltonians, considered 
suitable for describing mixed valence systems, are explored 
theoretically in this thesis* 

The introductory chapter of the thesis describes briefly 
the experimental background and the existing theoretical 
approaches to the problem. In most solids, ions exist in a 
state of integral definite valence. However, in a number of 
metallic systems containing rare earth ions, e.g., Ce, Ce AI2, Yb 
AI3, collapsed Sm S etc., the rare earth ion is in an unusual 
state of mixed or intermediate valence, % , such that 
n< n v < n + 1 where n is an integer. The mixed valence 
systems have many unusual properties, viz., absence of 
magnetic ordering (even though the ionic state with valence 
n may be magnetic) , anomalous resistivity, large specific 
heat, complicated (nonmagnetic) phases. The mixed valence 
phenomenon, specially marked in the early, middle and late 
lanthanides, occurs also in dilute solutions of these in 
metallic hosts, and in actinides, as well as some transi- 
tion metals. 



Theoretical attempts at understanding the phenomenon 
concentrate separately on two different aspects. First there 
Is the question of energetics; i.e., for this interacting 
electron’ system,, what is the ground sta+e ? Why is It one of 
nonintegral valence ? Is the transition (from insulator to 
metal, and from Integral to non integral valence) continuous 
or discontinuous ? Since correlation effects are very impor- 
tant, this serious electronic structure problem has been 
tackled so far largely through phenomenological approaches. 
Second (and this is the area of our concern) , there Is the 
problem of understanding the detailed magnetic, thermal and 
electrical properties of mixed valence systems. Here, the most 
common approach uses the microscopic model Hamiltonian intro- 
duced by Anderson to describe a magnetic impurity in a metal. 
The Hamiltonian is 


H = \ £ d + l 

o 


kcr 


+ 

kd a kcr a dcr 


h. c . +■ J e 
kcr 


U 


n - + l 2 n dtr n d- 

a 


k Hcc’ i 2 


da d 

where the atomic ^-electron (the subscript d refers to 
electrons in the localized level of the impurity — f-orbital 
■in case of ra re-earths or actinides) has energy with 

respect to the Fermi energy, mixes with amplitude with the 

conduction electron state of energy e k , and repels a d- 
electron of opposite spin with energy U . In a concentrated 
system, e.g., a lattice, there Is a d electron orbital at each 
lattice site. As e d (a phenomenological parameter) tends to 
zero, the states with valence 1 and valence 0 ( 'n|j> = 1 and 
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(n y - 0) become degenerate , and so the possibilities of two 
different valences and intermediate valence are contained in 
the model* Th^ model, however, dO'-'S not Include degeneracy of 
d or f orbitals, crystal field effects, etc* But it is the 
simplest one exhibiting two valence states. Earlier work on 


the model has concentrated on the r good' magnetic impurity, 


f» 11 l V kdl 2 P £p = A C p £;b , 

states at Fermi level) and 


= conduction d~ctron density of 
; ^ (symm^+rlc Anderson model) 


or |e d | >> a (Konds impurity). This and earlier Hairree* 


Fock work are briefly reviewed in Chapter I. 


start our work on the mixed valence problem with a 
thorough analysis in the Hartreo Fock approximation (HFA) * 
For a single impurity this was done by Anderson, and an 
important extension to include a d electron conduction-hole 


attraction (a Falicov Kimball term, Hr 


W a a do kW' a to'W has 


been made by Haldane. Haldan« found that when the parameter 
2 

C ( n» ¥ p ) is sufficiently large, there Is a region in 

e d U 

the parameter space ( , -£ space) where a nonmagnetic and 

a doubly degenerate magnetic states are simultaneously stable. 
In Chapter II w c investigate this region :»n detail, especially 
around * 0 , locate th« self consistent solutions to the 

H-F equations by computer calculations, find their energies and 
the coexistance region. We also obtain approximate analytical 
forms for the latter and show them to be accurate. We find 


basically that for C > % h there exist in the parameter 
space around =0 and = - U , two ’mixed-valence 1 
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regimes of width ^ C. Within each such mixed-valence regime 
there is a line on which the magnetic and nonmagnetic solutions 
have equal energy, but <n$ > and < s d > ( s d = na + - n$ + ) for 

the two states are very different; thus this is a line of first 
order phase transition. 

We next consider two °uch impurities interacting through 

a hopping term of the type X + ti-c • )» 3n HFA > 

oo 

the aim being to obtain the energy of interaction E} for 

oo 

various single impurity HP configurations CE^nt = Eg - 2E^) ' 
for the case C / 0. Denoting the nonmagnetic state as o and 

the magnetic state as a , we calculate the energies Ej n t $ 

00 00 0-0 

E^ n ^. , Ejf , Ej ^ where the superscripts, oo for example 
mean that cne of the impurities is in the state o and the ether 
in the state o • In addition to th«se diagonal terms there 
is the possibility of transition between individual ion H F 
product states, i.e., processes of the type ( oo > |ooi>. 

The amplitude for this is \ . 

The above serves as a preliminary for investigation of 

the mixed valence lattice (Chapter III) . The existence of 

three states, two magnetic and one nonmagnetic, at each site, 

is modelled by r.he eigenstates of an S=1 spin, with a single 

\2 

ion anisotropy term AS^ simulating the energy difference 

i 1 

between the magnetic (S z = + 1) and nonmagnetic (S z = 0) 

states. The interactions obtained above can be written as 
spin-spin interactions. The important points are that (i) these 
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interact! ons arise from states cut-side the restricted three 
level manifold-, so they should be included 3 r any complete 
theory working explicitly only wi^h these thr«e states and 
(ii) the spin - spin interaction* 5 ar^ obtained as functions 


of a d , C, U and A . We discuss the possible phases of this 


interacting ’spin* system, which has '-he mod'-l Hamiltonian 

H = i |S l2 + f J S^S l + t Xs|g^ + r iis 2 2s 2 2; - » 
i id 4 i J f i j 

>z > 0 a nd < S z > f . u i- ) 2attice 


The orderings <S%> ^ 0 and < S *> r ; iM fl++n - £ 0 correspond to 


ferro- and antiferromsgneti c phases'. We obtain parameter space 
conditions for this to occu^, and mean field estimates of the 
transition temperature T c . In addition there is the possibi- 
lity of < s| > / 0, which is a nonmagnetic coherent valence 

fluctuation phase. One important result we get is ‘■•hat this 
new type of phase is more stable than the magnetic phases for 
a sizeable range of values close to zero. This is in 

accord with experiment. We compare transition temperatures and 
briefly discuss the nature of ■’"his phase-. 


The HFA is correct in the gross, bu% for a single 
impurity for instance, misses the Hondo effect-, which is 'due 
to transitions between the HF states owing 1 o fluctuations 
beyond the mean field,. Due cc the Anderson orrhogonality 
catastrophe, such transitions are nearly blocked, r.he Kendo 
temperature Is small, and the HFA good except for T < T K . 

A formalism taking precisely rhe above fluctuations into 
aocount, for the symmetric Anderson model-, was discusspd by 
Hamann . We extend (in Chapter IV) this path ‘integral formalism 
to the nonsymmetrical Anderson model, and show that its 
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p rtit^on function is the samp as that of the Kor-do model 
as, obtain^ v 

Dy Yuval and Anderson. - From this analogy w= can 

calcu late* t 

eff &s a function of . In contrast to the 

Schneffer^Woiff va i ue J = ” (—■- “hq;) > do ^ s net 

dlv^Tsp a«= u u 

' * e <i ** 0 , though it dees increase to a large value* 

ThUv, the ssymmetyi c And'-rsoo mod^l tends to a strong coupling 

Kondo system with the Kondo temperature T^. -*■ a as 0 * 

W,. f , tally d i g cus s the impurity with C ^ 0 in this formalism, 

and s^ow that the partition function is identical with that 

Of a t:Wn T»^-t 

■ Pb °nant. level model. The la ■'‘■ter, unfortunately, is 

not soluble n , , 

^ Un contrast t.o ^he one resonant level model 

the Toulon 

mu 0 

ft^thod of approach of -ihe previous pare graph cannot 
easily be carried through for mor^ than on-' impurity, e.g., 
ttice of mixed valence ions. The mixed valence problem is 
characterized by several energies,- e.g.,- "the Ko-idc remperat-ure, 
the energy djff Pr(?ncn b^twe^n two diff®rsE ir valence states, 
the energy of Interaction between ions, be’ag nearly the same. 
This limits the utility of HFA in which T k = 0* The path 
integral method is tractable only for E-j n j- = 0* Thus we need 
a mo .hod which includes all ■•'he so effects. Preliminary results 
using such a method are presented in Chapter V. We use Bloch 

«md De Domini cis formulation of expressing the partition 

function 2 nf a , A . . . 

°r an interact •. ng system as a product 


where z 


z = z • z 
o Vp< 


'q.p. is, la appearance, the partition function of a 



"quasi particle" with discrete real energy levels. For the 
Anderson model with U = » , we find the energies of these 

quasi particles (three in number, one nonmagnetic and two 
nonmagnetic) within a reasonable approximation, their dependence 
on temperature and model parameters, and h^nce thermodynamic 
quantities, viz., the static magnetic susceptibility for all 


temperatures and , the valence fraction 

» 


(defined < n^> 
to be 0( ) , 

ground state at 
fraction of 1 - 


). The zero temperature -susceptibility is found 
which Is consistent with other estimates. The 

zero temperature is nonmagnetic with a valence 
T 

•J& showing that the ground state is a mixture 

A 


of the two magnetic and the nonmagnetic states of the non- 


interacting system (the term 




c do + 


having been 


taken as interactl on) , with a predominant proportion of the 
former. This corroborates the singlet ground sts+'e idea and 
quantitatively matches the estimate of Yafet and Verma. 


It is possible to Include the C ^ 0 and finite U cases 
within the formalism although we have not gene into the details 
of these. For the two impurity problem also it is possible to 
find out the energies of interaction betwr-n quaslparticle 
product states, and tackle the lattice problem via a pseudospin 
model. This pseudospin model is, however, very different from 
that obtained from the HFA, as the magnetic and nonmagnetic 
states themselves are different. The problem Is treating two 
impurities in this formalism arises in proving the fundamental 
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relation Eg = 2E^ (for E in ^ = 0, or two nonin^eracting 
Impurities) . 

Lastly, we must mention That the outcome of The above 
treatments are not yet in a state to be compared with experi- 
mental data. The theory or the models, as already mentioned, 
are too simplistic and therefore comparison wi^h experiments 
would be premature. Bu* we expect these models and their 
treatments to provide a basis for further improvements and 
make the theory more r^&lis^ic. 



CHAPTER I 


INTRODUCTION 

Seme rare earth compounds such as chalcogenid es , in^e metallic 
compounds, Cerium metal, and more recently, some actinides have been 
the subject of a vide range cf experimental investigations for sane 
time. The Interest in rare e&rt.h (RE) chalcogmdes such as SmS, 

SmSe, TmSe, Tm Te etc., arises from the fact that although under 
normal pressure they are semi conducting and usually in a 2 + (atomic 
configuration 4f n ) valence state, on application of hydrostatic 
pressure, they undergo an abrupt is ©structure 1 change with a large 
change In the lattice constant - , , to a conducting ste*n with an 
unusual valence In between 2 + end 3 + (configurations 4f n " 1 5d x 6S 1 "’ X ) 
state. In this phase (which is also called a "collapsed" phase) 
lattice constant corresponds to the weighted average cf those 
corresponding to the 2 + and 3 + states. It is possible to find the 
ratios of fractions of 2 + and 3 + states quite accurately because In 
lanthenide series, due to lanthenide contraction, lattice constants 
vary very systematically with the valence or number of electrons 
in the atom. This metallic phase (collapsed) shows anomalous 
magnetic properties — saturation of the moment at low temperature 
and absence of magnetic ordering; resistivity minimum; specific 
heat anomaly, and so on. The intormetallic compounds ce ALg, Ce AI3, 
Yb Alg etc., and the actinide metals such as Np, Pu and some 
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and Pu A1 show 

behaviour similar to this collapsed metallic or mixed valence 
(MV) phase. The mixed valence phenomenon is also seen .in numerous 
metallic hosts containing RE impur'Mes . 

The interest in irt^rme.li ate valence or fluctuating valence 
compounds has be-'n signifi care end is growing as Is evident from 
the large number of recent publications and conferences on the 
subject. Since extensive and upto date review articles'*' exist on 
both experimental and fch^oreti cal aspects of the subject we do 
not consider it necessary to attempt a detailed review here. 

We will, however , discuss the salient experimental features 
wherever necessary to explain or justify our choice of models. 

Theoretical attempts understanding th* mixed valence- 

phenomenon concentrate on two different aspects. First there is 
the question of energetics i.e., for this interacting electronic 
system, what is the ground state ? Why Is it of nor? integral 
valence ? Is the phase r ran i or continuous or dfs continuous ? 

In this electronic structure problem, correlation effects are 
very important as, clearly, the extremely localised f-orbitals 
are involved. This problem lias so fer be'~r J cackled largely 
through phenomenological approaches and we have not gone into 
these, although at ^ lines we have based our models on these 
approaches. The second aspect of the problem forms the area of 
our concern. It is that of understanding th: detailed magnetic, 
thermal and elect rical properties of the mixed valence system. 


compounds of U, Pu such as USn 3 , UAlg, WpRhg 
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The magnetic susceptibility , specific heat and resistivity 
of MV compounds are similar '‘'c rho*? of a Kha^c system. This, 
and the lack of magnetic ordering, indic&t* '"■hat- '-nt^rsite inter- 
actions are extremely weak and the study of a single MV ion will 
take us a long way towards understanding rh-- whole system. 

The physical properties which we seek •< c expls-”' are Inter- 
mediate valence (f luctuatl on or mixing), static magnetic suscep- 
tibility and magnetic ordering. There has b^en some speculation 
in the literature about the time nature of mixel valence in 
these compounds — whether valence m a sample .is spatically 
homogeneous and fluctuating in time or whether different valence 
ions are spatially distributed like an alloy. Th^ dynamic probes 
for this are Mossbaur effect, neutron diffraction and XPS. The 
first of these methods gives oily an average valerce but the much 
faster probe, XPS, clearly shews the two valence states. Thus it 
is established beyond doubt that these are actually fluctuating 
valence systems with chars c* eristic time typically of '’’he order 
of 10” ^ seconds. The typical magnetic su s cep-” ■* bill *y versus 
temperature curves have beer. °how;n in fig. 1(a). They are 
cheractsrised. by a milt* peak b->low which the susceptibility is 
slowly varying and at z-'-vo T.smprro+uTr has a finite value (non- 
magnetic behaviour). A* much higher ■ emp~rature the susceptibility 
shows Curie -V/rlss like behaviour characteristic of a magnetic 
moment, 1 Thus, while the high temp-mature susceptibility does 
show th^ ion to carry magnetic moment, all search for ferromagne- 
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tic ordering has been futile; only one compound TmTe shows 
antiferromagnetic ordering at low temperature; and Ce AI 2 i-s 
thought to do so, though not- conclusively shown yet. A typical 
specific heat curve Is show n in fig. 1(b). Zero temperature beha- 
viour of these properties suggest heavy Permi liquid. 

Any theoretical model for such systems must- start wi*h a 
narrow f -bend (4-f for RE 5-f for actinides) or f-stete in case 
of a single ion, and a wide conduction band. Then the strong 
correlation between highly localized f-s+atp must, be taken into 
account. One such model is the Anderson model for localized 

p 

impurity in metals and we start with the simplest non-degenerate 
or single orbital Anderson model, the Hamiltonian for which is, 

H = t ^ °ka °ko + 1 £ da c d cr c da 

ka a 

+ J v kcr (c kr" c da + c da c k a ^ + 2 \ n da n da ^ 
kg a 

The first two terms describe the conduction band and the 
ionic local level (subscripted »d 1 for historical reasons but 
actually describes f-ststes) end the last cr.e correlation 
between local electrons. The third term describes a hybridisation 
between local and conduction states and arises from the fact 
that uiesuLved in a solid ? the local level is no longer ortho- 
gonal tc the band levels even if originally of a different 
symmetry. The single orbital case is not too unrealistic because : 
in systems under investigation, due to large Coulomb correlation j 
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(U) , cnly two ionic levels are close enough tc each other and to 
the Fermi level to matter. It however gives wrong megnotic moment 
values. But we believe that the essential statistics and quali- 
tative features will not be damaged due to this choice. In RE 
systems under consideration, the coulomb repulsion U -v 5 to 
If) eV, the mixing term V is responsible for the width of the 

g 

local state (A = ir p 6 | V | ) and typically we have A -v .1 eV. 

The distance e<g of the local level from the Fermi level 
decides the extent of valence mixing — closer it is to the Fermi 
level, more energetically favourable for It to empty into the 
conduction band. Therefore for valence fluctuation phenomena, 
the case of interest is when £ d is close to the Fermi ’level 
or l-e^l ^ a . The Anderson model is a strong candidate for MV 
systems also because the properties cf th-^se systems are dose 
to a Kondo impurity and sc are those cf the Anderson model as 
shown by S. W. transformation. The S.W. transformation, however, 
ices net apply to |e^| < a the region of .interest here. The 
correspondence to Kondo model, therefore, does net become clear 
right at the outset. It becomes clear, in this work, only in the 
Path Integral treatment (chapter III) of the Andersen model. 

But, as is well known, the treatment of the Kendo model in 
time representation (the Path Integral method (PIM) reduces the 
Anderson model to this) takes only spin fluctuations into account 
and valence fluctuations dc net explicitly enter the picture. 

This is because in the HartredFock approximation (HFA) , which is 
the starting point of PIM, ’"he two solutions of the Anderson 
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model are cf the same valence. Tc accomodate fluctuations 
between different valence spates in the HFA we incorporate into 
the Andersen model a term 

g C I c ka c k f c ) ^ ^dcr ••• C<2) 

k ’kcr a 

This represents a Coulomb repulsion of the conduction electrons 
by the localized electron. This is the in Fallcov- 

Kimball model 4 ' responsible for the phase * ran sit ion to a different 
valence state. The bracketed part in (2) can be represented In 
terms of Tomonaga bosons. (2) then represents coupling of the 
local electrons with a boson field. Thus coupling with bosons 
such as phonons can also be treated in the same way. 

In Chapter I we have taker up this modified Anderson model 
in the Hartree-Fock approximation at zero temperature. For a 
single impurity this was done by Anderson 2 who found that- there 

A 

exist two regimes in the parameter specs ( -p); one for which 

<n(jt> = <n^+> signifying a nonmagnetic regime, and another 
for which <n d+ >^ <n d+ > , signifying a magnetic regime.- The 
implications of this phase diagram and transition are discussed 
briefly in section 1 of chapter II. In the next section we add. 

Fall cov- Kimball like screening term to the mo^el as dene by 
Haldane 5 , and proceed with its Hartree-Fock(HF) analysis. We 
concentrate on tb/» large U, large C case for varying e a . We 
write down the HF equations, and their stability conditions, 
Haldane found that when the parameter C Is sufficiently large, 
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£ A I TT 

there are two regions in thf parameter space ( around 

6^ = 0 and 63 = - u, where a stable nonmagnetic solution co- 
exists with a stable magnetic one. We investigate these co- 
existence regions, especially around £3 •= 0, n detail by locating 
the solutions to HE equations by computer calculations. ¥e show 
that the coexistence region Is of width v^\ C. From the calculated 
energies and d-occupati on numbers ( ) of 

the solutions we show that within this MV regime (coexistence 
region) there is a line in the ( €2. , 2 ) space on which the 

A 

magnetic and nonmagnetic solutions have equal energy, bu* differ- 
ent • Thus this is a line of first order phase * ransition. 

This phase transition is, of cours-, spurious, although it is 
believed that the energies are of correct order of magnitude. This 
spurious phase •'tans It Lon disappears and gives rise to a gradual 
change from magnetic fc nonmagnetic behaviours when on" takes into 
account effects beyond the mean (or HF) field. This has been done 
In Chapter IV«- In section 2 of chapter II we find appropriate 
snatytical expressions for the energies end boundaries (in the 
parameter space) of magnetic and nonmagnetic solutions for large 
C ( 7T£ ) and show them to agree quite well with the computer 

calculations. 

In the next section (sec. 3) of this chapter, we consider 
two impurities each described by a modified Anderson model, and 
interacting through a hopping term of rh° type C^g + h.c., 

as done by Alexander and Anderson. 6 We find within th* HFA at 
zero temperature the energies of interaction where = 
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Eg - 2E-j_, for various configurations of the two impurity system. 

Here Eg is the energy of the interacting ions, E^ is the HP 

energy of a single non! rr eract-t og ion, a nd by " configuration" , 

e.g., denoted by superscripts o a in E° ° , we mean that one 

of the ions is in th^ nonmagnetic state en A *he other In the 

00 

magnetic state c • Thus we find the interact! m energies Ej_ n ^ , 
E int > s int anc ' E int * In to these diagonal terms, there 

is the possibility of moment transfer or exchange as a result of 
interaction. We show that the processes of the *yp« |0cr>-*- | a0> 
have an amplitude x • The results of this section .are used 
in the next chapter to construct a Hamiltonian for a lattice of 
such interacting ions. 

Thus, in Chapter III, we consider a laotic^ of mixe^ valence 
Impurities, the preliminary energetics of which are provided by 
the HF treatment of the previous chapter. In chapter II we show-- 
that each impurity (with C ^ 0) in the MV regime can exist in 
either of the two degenerate magnetic or the nonmagnetic HF 
state, the magnetic an^ the nonmagnetic stages being of different 
energies (generally). This situation is simulated by representing 
each ion by an S = 1 spin and adding a single Ion anisotropy 
term, Asf , where A is the difference between the energies 
of the two HF states. Now the interactions obtained in Chapter II 
can be written as spin-spin interactions. These interactions, 
say between an ion at site i and another at J , are cf the type 

J ij s z S z » B ij S z s z » x ij s x s x 9 ••• » 

where all the parameters can be obtained from the results of the 

last section of Chapter II. Thus we have the Hamiltonian of the 



MV lattice as 
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1 i 3 v i 2 3 2 r 
^ = I Ji jS z S 2 + I BijS z S z + A I 

13 13 


i 2 


s z + * l sj S J + M 
1 ij * Z 


( 3 ) 


Next, in section 2 cf Chapter III, we consider the various 
possible orderings of the above Hamiltonian within the mean field 


approximation (MFA). We fin'" 1 parameter space conditions and 
transition temperatures T c for the various orderings, ferromagne- 
tic ( <S Z > = 0), antiferromagnetic ( <S Z > / 0) and a turd 

sub lattice ■ 

kind cf ordering characterized by <S Z > ^0, <S X >^ 0. This 
last is the coherent valence fluctuation (VF) phase. On comparison 
of transition temperatures w« find that this nonmagnetic VF 
phase dominates over most cf the temperature range upto very low 
temperatures, it very low temperatures, w^ .expect Kendo fluctua- 
tions tc t&ke over which makes the ions nonmagnetic anyway. 

Thus it seems to bear cut the experimental results cf absence of 


magnetic ordering in MV systems. 


As has been already pointed cut, the HFA is co-rrct in the 
gross, but, for a single impurity for instance, misses the Kondo 
effect which has beer, shown ~o exist for an Anderson model. The 
HF wave functions which are eigenfunctions of different 
effective one electron Hamiltonians, are not orthogonal to each 
other though nearly so as shown by Anderson 7 orthogonality 
catastrophe. The small overlap b-'+woer these states causes 
fluctuations, which are beyond the mean field approach, and 
gives rise to Kondo effect and a vanishing moment at very low 
temperatures. To take these effects beyond the HFA into account t 
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O 

we follow Hacann's path integral method which ho user for the 
symmetric Anderson model. We extend his method tc the case cf 
asymmetric Anderson model (with C = 0) In Chapter 17. The method 
consists of choosing hopping paths between the two HP solutions 
which correspond to minima ir the free energy functional, and 
calculating their contribution to the partition function. We show 
that the partition function Z has the same form as ^hat obtained 
by Yuval and Anderson 9 for the Kendo model. From this analogy we 
find the parameters J, % of. the equivalent Kerr's model. By this 
method wg are able reach upto = 0 which the Schrieffer- 
Wolff transformation fails to ^ as it breaks ^ own as | | < a • 

We shew that the effective exchange coupling J e ff does increase 
tc a largo value but remains finite as |e^| 0 , taking the 

system to strong coupling Kendo limit-. In section 4 of Chapter IV 
we turn to the C ^ 0 case and shew that the Z for this, when 
scaled following Anderson and Yuval^s method, reduces to that 
for a two resofli£u + ' level model in contrast to the one resonant 
level model obtained by scali ng of the C - 0 Anderson model# 
Unfortunately, the two resonant level model is not soluble and 
therefore this method cannot give physical quantities for the 
C / 0 case. 

This path Integral meshed cf Chapter IV cannot be easily 
carried through to describe a many impurity situation either; 

The treatment of a la-*-ice of MV impurities in Chapter III is 
meaningful only when Kendo fluctuations ere co a much smaller 
scale (Tk *v» o) compared tc I on- ion interactions# Aroun^ = 0, 
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these various energies, the Kondo tp 2 p-.rat-ur= , energy difference 
between different, valence scenes, energy of interaction between 
ions, are all of the sam* magnitude. This lfmit-s the utility of 
the HF approach chapter III. These difficulties prompt us to 
investigate another method which holds ‘"'he premise of accomodat- 
ing nonzero C as w^il as many impurities, for all parameter 
ranges. Chapter V presents preliminary results us-o.g such a 
method, the quasiparticle method. 

In section 2 of Chapter V we very briefly i <- scribe : 'he 
formulation cf Bloch and De Domini cis, of expressing 3 of a 
large Interacting system as "hat of a * quasipar"? cle 1 with 
discrete real energy levels. The self consistent equations 
determining the quasiparticle energies for a U = « Anderson 
model (C = 0) ere approximated upto certain terms in section 2 
and the Implications of these a pprcxtmati *>ns, discussed . 1 itoong 
the three quasi parti cles , two magnetic and one nonmagnetic, we 
show the nonmagnetic one to be the ground state- i- + zero and low 
temperatures, while at higher temperatures th° magns^lc states 
lie. lower. Once the energies of these quasi particles and hence Z 
has been found, one can calculate thermodynamic quantities as 
derivatives of Z , in sections 3 and 4 wo fv*,d tfco static 
magnetic susceptibility, and valence fraction 
C < > ) . Numerical values have beer obtsrnei by computer 

calculations. The results on susceptibility are compared 
with the renormalization group calculations of Krishnamurti, 

Wison and Wilkins 10 , and it is found that around I e I < a 
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our magnetic susceptibility curve lias tiie come general shape 
as their calculated ones, although our results are higher by 
a factor of -4 (not uniform over the whole range of tempera- 
ture). This may be uue to the approximations involved, or the 
difference between a U = D/2 (Kin!) and u = 00 (this 'work) model, 
lor the s , < 0 , Is, I » A case we show that i\ matches standard 
results at zero and high temperatures. So we conclude that 
the approximations involved are good. 

In -.he lest, sechl x, 5 , of -irs chapter , w- briefly discuss 
further pcs si bi lie ‘cs ^ o f-h' method® The C / 0 c&sr ^o<-s not 
Sf=m in h-e enable . ' Th^ i betw-’t •■w: impurities is also 

possible to *p&1 with. The difficulty in this arises because of 
the formal difficulty ?.n shewing Eg = 2 E^ for c pai^ of n"n- 
1 n t e r & c t i r - g lens® Bu^ W ' 3 b^liev" "hat ah' ir.'~ * no ■"* is cu pab l r of 
dealing wi*h these problems ac* h.-ic« th«.+ 'f a. le-'-^i ce of MV 
impurities. 

Ir *ho concluding chap' •t w- summari. z? tar results Obtained 
in, en H wisdem gained fr~in chapters II thrr ugh V. P- ’Slly, w ? 
suggest ■■hat. although 3 ** *hr p^snir ? T fr r f , results can only 
be compare ■ quf ll.ta^i v^ly wi-h exp^'m^n^al de^mcre realistic 
features like multiple c^bitals, crystal fi^l^e etc®, c&r. be 
incorporated "ito rhr mo ■’•'l 1 " ■"'■^■" cement ‘~f chapter V* then the 
results on physical quag- --s will br mev tpp for comparison 

wi^h experimental data. 



CHAPTER II 


THE MIXED VALENCE PROBLEM IN THE H*-RTRES-FOCK 


APPROXIMATION 


In s chapter we cpply tie simplest, th~ Ha ’"tree. Fcck 
approx’’ me ei on (H FA), to t h~ -tenders on model and. tt c 
vptsj o'ivccc.-'-'S -by Haldane 5 C G 0 ct-M?), £"? ' ~h?u a pair 

of such !:is, 

The HFL -n suits for the Anderson model fare well known an- -1 
''a section 1 w- b-iefly describe. lr. 

In sectfon. 2, w<- tuk~ up th r mc^lfied Ahdensor model 
(G ^ 0). l, h'th HaHan^'s w^k as background, we discuss In r '^^s.f 1 
+ hn mlx^s v&l-mc^ regim ! ~ in the per&mr "er sp&c 3 . W-> flnt ^h* 3 
bourse rl ►•s cf “Ills reg-m°, ^he er-rgi'.^s rf th^ two solutions, 
art 1 th^-ir v&l^nces by computer calcula tdons . Then, w~ find these 
quart - ! •“ f-s for large C (G > A ) in app—xlraa + ••> arely-ical form 
and shew th^m to be close to numerical result;--. „ 

In srcti:n 3 of :h's chapter, w<= consider*, in the HFA, 
two such Impunities interacting via a hopping trrm -f ~ho *ype 
Xc ri^ Q d 2 + 11 • c ' Farms 11 has pa really address^# himself t r this 
work. W~ analyse chr problem following Alexander an'' 1 Anderson's® 
HF approach and find ohe various en^rg* cf inhere ctj on. 
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S»c«*f 1 : THE HNDERSON MODEL IN THE H.-RTREE FOCK .aPPR QXIMHT I ON 


The- nc.‘ti' • g^ r ^ra^s lo ■* r- r s r. i mo*' d has b; 
HFA by ilrh-Tsoc h'ms^lf « Tli~ Hairh l^or* &~ 


n the 


, y 


H ~ J £ lze c k a°lza + \ c " 

1 CC 7 O 

+ (, , 

ucr key 

op-res.- :.*is cor dc 'i ~i~ -l- ctr~“ , s * . c w •' s' ben*', wv*h ■v-g^s 8^ 


6' l :: 'cLa 11 d--o-' i '^ / kd' C kcy c da * 
° + c" r Cl ) 

t • re 1 <"rr y 


mixing wi'tnh ""he lccdiz" dipu^' oy level "f e -'"gy (•-•rrrrg-: =s 

measure* from +h-^ F°nrH l~vel) vs a nrxing ''* r v® -f g-r-^gs-v, y^ 
whi l-"- local *'lo errors of epposi x? spin r^p-l --ech '■"h" ^ win 
coulomb on orgy U> I* th~ HFA, each local ?loc*r^r of spd cr 
:ts feoppos«r1 be b* governed by rho Ham' 1- "nt an 

I *k c k0 C H0 * (E d + I <J V. 0 >)n da * j ( c ka c d0 •*• h - c - ) 


(V is an av-'-rag" m*x 4 *g strength, (V k ^) av )« 

This ?«snen* level Hear l v .e it an has well k~„cwo 3clu-,t~rs and 
gives, si zero ^emp-re+u^-', 


<n do> = / P d ff (e) <1e 


, p 




where p do ( £ ) = £ ~ t 

( 0 *■• £ .. -- — < 11 , >) 
v.iO Z CL- >0 

with A = it£. 1 7 1" 

Collect *r.g exprosei •'o.s for local ©l^c-r^rs ;f bo^h spn~? one 
obtains - v he edf cons* s~~z ,i ’ '-que^i.ins , 


< rs. 


’♦ > = I 


t -x e d + CT<n dr 


< n . . > = A. cot "' 1 


-!- U <n, , > 

U Cl+ 
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This set cf equators he s no types :f s olu - J' ere , magnate, 
with < > 5/ > and nr “me gr •' 5 c wh=n <xt, + >- =. -< rj <a+ > « 

The magnet c solution* exl s fr rely wh-^n U i ^ • When this 
ccn' 5, i* ; i t. * s sa*-i sf 4 * r, -he Eicgn~- r ?c ecluf. exi s-s only within 
a rang-~ of a^ sp-'-cifi by — U £ £3 < 0 • WfMr ~h?s range, 
the ncran&g , ne > t c ft clu -”3 t 1 s tnsbahl^ (stability of s ilu- f ons will 
he ^‘scussH ri -he Hex- s^c 1 *?.::) , rnd cu*si only the non—. * 
magnetic one exists,. Thus w-* have he phas -■ diagram when in 
fig.* 2* 

This phase dj a gram shews a sharp trarsi" 1 or b : -w -•-•• 0 
magnetic an-* nonmagr-ttc chart at p&r "leu 1 st pcr~ - e ir the 
phase space* This, however, should no* b~ ^sk°;. lit -’'ally. K^ep- 
■*ng in view the character cf HF-s., which pr^clu^s slow oscilla- 
tions between, say, the -wo tego'-'era-e megn-^ic ? eluti trs. .such 
oscillations are ch~ heft of ’•hr Kondo offset ir. A er** near iy 

V 

blocked by -ana ersc-n f s cTthogcna 1:? -y catastrophe according ?c 
which the matrix element be-w^o tw c detetmin^n-tal wave f uoc "i ns 
cons- 1 ' meted from ff otti*' scattering s ,w £tes .(as in the HFi. for , 

An 'i ere or model for local ol-cir - 0 of ~ach spin, for a magnetic 
solution) is- vanish! igiy small when the phase shift? at the 
Fermi surface of the tw r scatter*;? ng states arc greenly , 

d?ff ^r^rrt*, Thus it * he HFA we miss ! "he Kondo off ^c" . However,. 

1 *• gives us ce-r-ct. energetics s~.- p within th° magnetic .~r non- . 
magastic states Cf&r away from the phase trans? it or lln^), in 
r hat 'ho -ordinary theory cf .norma grot* c- - 1 mpari*- 1 •?> applicable 

in the nonmagnetic regime and in the magnetic regime-, Keroo 



18 



Fig • 2 Phase diagram for C = 0 Anderson model. 

(from Verma ) 
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effect Is regc. tiled If upci the HF s clip- ions we Impose the above 
m.-nt.t Tied oeci llehicns. Hrice the ’mpo^t’tce of <.r. HF analysis. 
These pof nos have be”" discussed o rietsjl by 3 Inn din. 


Section 2 ; ANALYSIS OF THE G ^ 0 ANDERSON HAMILTONIAN IN THE 
HaRTREE-F OCK aPPRCKIMaT I ON 

4 

The Ann rs ■ -> model has br^r. onc'-rily g^r - ' tali zed by Haldane* - 

to irtclud ” local electron conduct! ci electro t c^ulmb repulsion 

by l^clu^irtg in *ho And— son Hamilicr-isn s *erm g C I na<y) 

cr 

C J, , d kcx c k ! a' )• I " h£s b ?e". shew- ihat by *tk1ng terms upto 

lx y Xi. j/ O’ 

sec Oft order in g (the approximation with;” which Halsr? find 
Ri Sf-borcugh 13 have worked), an effective Hami l'-o^i an can be 
written as follows : 


: 8 da t d<sr c da + 

l c L°ka + 

U n , . n . , 

at at 

£ ( -n -i- > 2 

+ V I ( c, 

i _ k(J 

key 

, H- 

c -l- c . 

a a cia 

C = 4 p 

P 

« 2 



. CD 


The above type of y-' ,, &ch " - r. also appears in th° case of 


? rf er&cti en of th« local ”l-c *t-"s will s Boson field. 


13 


Wo now analyse the Hamilton an (1) In the HFA. 

Section 2(a) s THE HF SELF CONSISTENT EQUATIONS AND STABILITY 
CONDITIONS OF THE SOLUTIONS. 

The effective fl -In’ ss^r- by a single d electron of 
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spin a is, 

x 0 = edu+Ctf-C) G <n d0 ' > 

Thus for the two spins w* wri^-o 

x -f = 8 d + CG - G) < n_ > - C <n + > ... ( 2 ) 

x - = + (IT - C) < n + > - G <n„> ... (3) 

The average occupat.1 ~rs <n + > arc derived from the effective. d— 
electron Green's function as 

. _ A r° dw 

'V = % J A cirr^ jrrTs 


A 

TU 





... C 4 ) 


Equations (2) and (3) with < n + > given by (4) constitute the HF 
equations for the two spins, to be solved self consistently. For 
C * they reduce to the S'* found for the Angers on m Vol (sec- 
tior 1 ) % 


In this case also two types of sclu-ions exist - the non- 
magnetic ones, for which x + = x_ and the magnetic ones for which 
x+yb X f Th= conditions and regions of parameter space in which 
either type cf solution exists, or both kinds coexist, will be 
investigated in de+all m section 2(b). Before that w- lock into 
the conditions for s+ablll^y cf i-hrs- s olu^i '"'•3. 


The self consistent equations (2) and. ( 3 ) do net. indicate the 
stability of their solutions. For this we must first write the 
free energy in HF.a j- 



cud. cu 


(«- x + ) 2 +b 2 


A 

% 


I 

—A* 


o 


.0)* 


weld) 

X ) ^ + 



(<n + > + <n 


>) 


2 


U <n, ><n > 

-r -* 
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The ti -ps*-. two terms represent er-Tgios cf *he two resonant 
d-el~ctmns and the lest rw*' these :V hfc* are left; af-'mr part of 
the correlate on terms have been absorbed ;n -he “ffect'va d- 
el-ctron Hamiltonians. The V-ttvad"i m of Fjjp and the matrix of 
its derivatives te ^escribe- 5 4 r appendix 1. I 4 , *s clear from (5) 
that Fgp is a function cf the variables x + and x_ . The 
extrema represent the HF solutions. The -tabii.imy conditions, 
viz« t the conditions chat rfcrse ar? minima, are derived in the 
appendix and are. 


2 + C { nt(x + ) + n f Cx_) } > 0 • (6) 

1 4- C { n»(x + ) + n»(x_)l - {(U-G) 2 - G 2 } n»(x + ) n‘(xj> 0 ...(7) 

whelms n*(x. ) = 1 C 2L - tar ." 1 ... C8) 

n 2 & 

It car be easily shewn , by pu*! f-*g 4 -h" LHS of (6) equal to 
zero into (7) “ha* now (7) is nqr sc*i «f • sd , i . , (6) does not. 
imply (*>1 G" the cth^r hand, (7) implies (6). T hoofers, (7) is 
a more severe cr.ndi hi on and her- eft -m wo shall c ■ is > ' ~r only (7) 
for axami-il rg *h<- stability of the solu^i ms (assuming all the 
time | U > 2C) . 


Section 2<b) j FINDING SOLUTIONS TO THE HF EQUATIONS 

The HF equations are, fra equals’ 'ms (2), C3) » (4) &n>3 (8), 
x+. = c et + (u - G) n(x_) - C n(x 4 .) ... (2*) 

x_ = + CU - C) nCx*.) - C n(x_) ... (3») 


Subs tract (2») from (3 f ). Then 

x_ - x+ = U {n(x + ) - n(x_)} 

X + U n(x_) = x+ + U nCx*.) = K say 

NM» 1 


C 9 ) 


or 


♦ • * 
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We recall fran (8) that. n(Xf) are functions of the form 
n(x+.) = \ ( JL - ^an" 1 *fc) . 

So, graphically , solutions '*o equation (9), for either x+. or x„, 

z - K 

lie on the intersection of n( z ) and th* straight line - - } 
shewn in fig* 3* K is the intercept on z axis an' 3 varies from 
solution to solution according to . 

The solutions arc magnetic or nonmagnetic accor-’ 41 .eg as x+. 
and x_ are different or same. Ie the straight line - — , when 

K is too small (e.g., on the loft, lire Nq) it cuts the n( z ) 
curve at only one point, and therefore only on-", a nonmagnetic, 
solution exists. As K becomes larger (the line moves -,c the 
right) the line reaches a position when >t xs tangent to n(z) 

(dented line 1). This is when the magnetic solution becomes just 
about possible. We call it “he extreme magnetic or extreme non- 
magnetic solution. As the straight l"ue moves further tc the 
right, three magnetic solutions become possible (e.g., *hc line N3) , 
for which (x + , x_) = (x^, xg), (X2, Xg) or (x 1? x 3 )$ and also the 
three nonmagnetic solutions (x^, x^) , (x 2 , x 2 ) and (xg, xg) • On 
moving K further to the right, we again reach an extreme magnetic 
solution (dotted line 2), after which, again, only nonmagnetic 
solutions exist (e.g., l*ne N 2 ) . 

Por nonmagnetic solutions of the kind given by line Ni, 
i.e., when the d level is deep beneath the Fermi level, both 
spins occupy the level and th® d-el^ctrcr occupation number 
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is ^ 2. For good magnetic solutions (x+ and x_ well separated) 
the d- occupation is 2l 1* The approximate solutions x+ for the 
three regions (provided they do not come too close to the limit- 
ing ones) are the following : 

In nonmagnetic regime 1, i.e.., £(3 + U - 2C * A 

x + = x_ + U - 2 C 

In nonmagnetic regime 2, i.e., e - > 0 , x+ = x_ 

In good magnetic regime, x + ^ £ d _ c , x_ 2i £ d - G + u 

- .+ 

The range of for which above regimes exist, are dis- 

cussed in the next sub-section. 

Section 2c ; BOUNDARIES OF MAGNETIC AND NONMAGNETIC REGIMES 

In this sub-section, we use the stability criterion (7) to 
sift, out of the various possible solutions described In the 
previous subsection, those which are stable. This enables us to 
discuss stable magnetic and nonmagnetic solutions, and their 
possible coexistence fo ^ certain range of values of £ d * G e ^ c * 

Boundaries of Nonmagnetic Sol uti ons : 

For nonmagnetic solutions, x + = x_ = x (say). The condition 
(7) becomes, for stable nonmagnetic solutions to exist, 

1+2 Cn'W - U (U - 20) {n’(x)} 2 > 0 

or, { (U - 2C) + n»(x) U (U - 2C) } { U - n f (x) U(U - 2C) > > 0 

... cio) 

We have taken U > 2G and n T (x) < o. So the. second factor in (1€) 
is > 0. So we must also have the other factor > 0, or 
1 + U n*(x) > 0 or | n «(x)| < 1 , 
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for stable solutions. Clearly, this holds for all nonmagnetic 
solutions tc the left of extreme solution 1 (sgp fig. 3) or to 
the right of extrema solution 2, i.e., for £^< (>) that ccrres 
ponding to extreme solution 1 (2). These solutions are 


n»Cx) 



A 

IX 


1 


X 2 + A 2 


or 


x ~ / U A “ A 2 



for U » -n; A 


Thus for nonmagnetic solutions to be stable, 


x > /au or x< - /au . With /au “ x > - / a if 

/ K / TX / 7X / 

nonmagnetic solutions exist, but are unstable. 


x 


The value of 


e a 


corresponding to the lower boundary 

as 


= - /aU ip given using e.g., (2’) or (3») 

/ ix 

e d = x» - (U - 2C) n(x) 


= - (— ) 


(5 + t an~l (AU)) 

TX '2 , 4 '"IX ' ' 

A , , ’"oAx , ~ /CA / m A > 


- U + 2C - 2C v ’~ + 0 ( A v 0 + 0 

Similarly, the upper nonmagnetic boundary, where x = +^/ JQg 
is obtained as 


: d = 2C,C + 0 (4,5-) 


+ o (8A , £4) 

u 'U vu J 


... ( 12 ) 


Comparing (11) and (12) we find that nonmagnetic solutions are 
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stable beyond the region bounded by, 

( £ d t (S„c+2 c f ) ... (13) 

Boundaries of Magnetic Solutions : 

It Is difficult to locate analytically boundaries of magnetic 
solutions. However we present here some arguments enabling us to 
do this approximately. We support these results by computer cal- 
culations . The computer results and conclusions drawn from this 
analysis are discussed later ir this section. 


Pig. 4 shows a plot of energy F^p (arbitrary base) of the HP 
solutions vs. , as obtained from computer calculations. Accord- 
ing to this, there is a region of £ d in which two magnetic solu- 
tions exist. We now show that the dotted portion of the line 
representing magnetic solutions corresponds to unstable ones. 

Por this it is sufficient to show that the condition for stability 


is the same as that for an extremum to occur in the curve, i.e., 


the condition —jf = 0 ==> ' minimum in F, or equation (7). This 

9s d 

has been done in appendix II where we first show that -rr— = r) ==> 
3e d dZ + 

— ay = 0 and that this condition is the same as (7) with an 
equality sign; the branch denoted by the solid line satisfies 


the inequality and therefore corresponds to stable solutions. 


Thus we show that there exist ranges of s d for given U and C 
within which not two but one magnetic, and one nonmagnetic solutions 
are stable. We call these the mixed valence (MV) regimes, since, 
as has been already mentioned, the valence or d-occupation 
numbers for the two kinds of solution are very different. The 
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extrema of nonmagnetic solutions have already been found. To 
find the limits of magnetic solutions, we must solve the equa- 
tion obtained from (7) 


1 + C { n'(x+) + n«(x_)} - {(U - C) - C } n f (x + ) n»(x„) = 0 


(14) 


We are able to solve (14) analytically only approximately 
and for large C. To solve it, substitute for n»(x + ) In (14) 
and we obtain 

ft a /_ TT/TTf';\»2/_2 

=0 ... (15) 


1 _ c &/7I _ C a/tc U(U-C) A 2 /-re 2 

x|+ A 2 xf + A 2 ( X 2 + A 2 ) ( x 2 + A 2 ) 


Suppose (x + , x_) is a ’good 1 magnetic solution, i.e., n(x+) * 1 , 

n(x_) 2l 0. Then we can neglect the last t c rm in (15) since x_»A 
and |x+. I * 0(U). In any magnetic solution forming the boundary, 
one of the levels has to be within + yAu/u- a 2 (the limiting solu- 
tion) • Suppose x+ is this level. Then the other level, x„, in 

a good magnetic solution must be such that |x_| » x+. Then (15) 
reduces to 

C A/-n; 


1 ~ p 2 

X+ + ^ 

This has the solution. 


0 


• * t 


( 16 ) 


X+ = + f™ ~ A 2 ' -v + i/LI for C » Air ... (17) 

To find out x_ for above x + , we use the relation (9), 
x+ - x- = - U n(x+) - n(x_)> • 

For 34 . = +V’~- this gives, 

\f( AC/tt) - x. = - U(V(Vrcc)-l + A/nx ) , ... ( 18 ) 
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after expanding the tan” 1 function involved In n(x) , to lowest 
order. The solution to (18) Is, 


x_ 




... (19) 


For x+. given by (17) and (19), n(x + ) and n(x_) are, 

n(x + ) 1 41 5 a(xj * 1 - A /[ U- ( u+c ) V ( A Ac ) ] 


and using the relation (3 f ) we find a $ • 


e <j = x_ - ( U - C) n(x + ) + C n(x_) 





- U + uf: 


TCC 


Co - 0) /A. _ 0 ( 


AC 


U 


* - U + C + 2 4" ••• (20) 

Thus (20) gives the lower magnetic boundary. The upper magnetic 

ijc 

boundary can be found likewise, taking x + = -V“~ from (17). 
The equation for the line In phase space bounding the magnetic 
region then becomes , 

e d- c + 2 *±(g-2 V™ ) for U, G » ttA ... (21) 


Comparing this with the expression obtained for boundaries 
of nonmagnetic regimes, eqn. (13), we find that the width (in e<j) 
of either MV regime is 

W MV 21 C - 2V~“ . -2 Cf^j ... (22) 

aquations (21) and (13) determine the boundaries of magnetic 
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and nonmagnetic stable solutions, within the approximations 
mentioned* 


section 2d S THE PHASE TRANSITION LINE 

As is clear from fig. 5, the HP energy curves cross each 
other within the MV regime. In the parameter space this repre- 
sents a Line across which the HF solution with lowest energy 
changes character, l.e., on one side it is the magnetic one and 
on the other it is the nonmagnetic one that has the lowest 
energy. Thus this is a line of phase transition* Around this 
line, the valence fluctuations are expect*/ to b--> the most 
prominent. Fig. 6, which shows a plot of valence or <n<j > 
against , shows that at, this line in parameter space, or 
the value of where the two solutions have equal energy, 

< n<j > for magnetic and nonmagnetic soluti ons are very different 
from each other. Sc this Indicates a first order phase transi- 
tion, unlike in the case of C = 0 Anderson model in which the 
solutions show continuous change in energy and < n^ > as the 
character changes and there Is no MV regime. 


To find the line of phase transit- 4 on we make linear 
approximations to the HF free energies of magnetic and non- 
magnetic solutions, given in appendix I. 


F HF = Un(x + ) n(x_J 


C {n(x + ) + n(x_)} 


2 


+ M 
+ 2 


In { 


2 

< + 



+ { n(x+.) + nCx„) } + 

+ const. 



-50 
F»S. 5 


- 20 0 10 £0 30 

— 

Energies of stable HF Solutions (arbitrary amount 
has been added to Fnp for each C) 
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Fig. 7 M.V. Regions in parameter space. 
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In the lower nonmagnetic regime, n(x + ) = n(x_) * 1. 

• • Fjjp 2L U — 4C + 2 £ ^ (23) 

neglecting the logarithmic terms. 

In the upper nonmagnetic regime, n(x + ) = n(x__) * 0. 

-^HF — 0 ••• (24) 

In the magnetic regime, n(x + ) ^ l, n (x_) ^ 0 

, » Fjjp — C ••• (25) 

The upper (lower) phase transion line is obtained by equating 
(25) with (24) ((23)) 

■So the equations for these two phase transition lines are 

found to he 

e <3 *” C + ^ = i ( ?r — ••• (26) 

Section 2e s SUMMARY OF RESULTS OBTAINED FROM HF ANALYSIS OF 
ANDERSON MODEL WITH C ^ 0 

Figures 5 through 7 show results of numerical calculation 
compared with our analytical estimates. We find our estimates to 
be correct, to within 0( A ) • For the display of graphs, the 
parameters are so chosen because, as mentioned in Chapter I, U 
is expected to be of the order of 'v 100 A • The parameter C 
is expected, to be much smaller, an order of magnitude smaller 
than U, £3 is a free parameter of f -he theory. 



The salient features found are : 

(i) Magnetic solutions exist only when U > irA 
(11) MV regime exists only when C > irA 
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(ill) There are two MV regimes, around 
and . n + _C + 2V~< E d <~tJ + SC 

•vC- 2/“ 

v n 


® £ £ d “v ^ 
wi+’h width 



of 


(iv) Energy differences be^w magnetic and nonmagnetic 
solutions in the MV regime are & 


(v) The phase transition within the MV regime Is a first 
order phase transition. 


Section 3 : TWO INTERACTING IMPURITIES, C £ 0 

Section 3a : THE MODEL AND THE SELF- CONSISTENT EQUATIONS IN HFA 


In this section we consider oho interaction betwooo two 
impurities, each of which is described by the generalized (C / £) 
Anderson Hamiltonian. The only ext-ra f-erm introduced is 
X° o, d2 a + h«n,. , representing the transfer of d-electrons 

between sites l and 2 due to wave function overlap. Another 
source of coupling between the impurities is via the conduction 
electrons with which the d-states mix. The Hamiltonian can 
thus be written as 


U 


H J 8 k c kcr c k0 + . ^ ^ s d c dia c dicr 
ka 1>cr ik* r . 

+ V L- (0 di<,°ka e ' + k.c.) + ( to d i 0 o d2 ff +h.c.) ...(ZO- 

0* y 1 j «m. 

where i = 1, 2. 


2 ^dicr n di“a^ 
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Wg Investigate this Hamiltonian in the HFA, i.e., we 

* 

reduce it to an effective one body Hamiltonian, 


H HF = l e k c lm c Tm + . £ ( £ d c dia c dio + L ' <n di-cf >n dicP 


+ 


k, a 


+V E (ct_. e 


1 s a 

ils-r. 


a , i , k 


'dia kcr 


+U.c.) + (X '= d l 0 c d2a- +ll - c - ) " li l <n dio ><n di-a > 

X ... (28) 

Thus the Hamiltonian reduces to that of two T>rscnant levels at 


sites 1 and 2 interacting with each other via th-" mixing term 

o 

proportional to x ...The quantities (<n^i a >) in (28) are to 
be determined self consistently. The self consistent equations 
for these are obtained as, 


n = (N, + N )/2 • • • (29) 

O “T — 

d 0 = ~(l/2V~ 0 )(N + - l\Tj { ( U-C ) d a - Cd" 0 > ... (30) 

where H + = 0/2 - tan -1 ((a d + (U-C)n" 0 --Cn 0 +V~ 0 )/A}]O 

c = + 5 n 0 = (<n dl0 >+ <n d2ff >)/2, d 0 =(<n dla >-<n d2ff >)/2 

V 0 = \ r [{ ( U--C) d "** 0 - Cd 0 } 2 +x 2 ] and |x 2 = X° + ^ T ••• ( 51 ) 

... ( 31 ) 

* d2 * s a transfer term arising out of the localized conduction 

electron mixing and Is given by 

ik* (r-^-r ) 

X 0 = -V 2 I ® — ... ( 32 ) 

12 k °k 


Diagrammatic&lly, it originates from propagators of the type 
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lc 0} 



1-tC a) 


20c 


V 


X 

y~ 


where dotted lines denote d-electrcn propagators and curly 
ones, those for conduction electrons. From the above diagram 
it is clear that we should actually have 


12 


(03) = V 2 l 


ik- (r-j -r 2 ) 


k “ " °k 

In the next order, i.e., 0(X^ 2 )°, the form (33) produces 
RKKZ exchange interaction via scatterings of the type 


(33) 


c ^2 


x. lc 
->k 


kc 


2c / 

V 


'v ej<s, 

/ ' s X 

lc k'-c 

/ 1-0 2-c 

c 

However, for the sake of simplicity, we omit w -dependence of ^j_2> 
use the form (32) and lump both kinds of hopping, direct and 
through conduction electrons, into X and hereafter use 
only X , The four equations (29) and (30) (two for each cr ) 
have been obtained in a straightforward manner following 
Alexander and Anderson^ and are the same as these found by 
Vanna*^. So we do not give the derivation here. . 
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The general solution of th^se equations is quite difficult, 
and we consider only the case of weak coupling between the 
impurities, i . e . , retain only the loading order ^orms in A • 

We further assume that the solution tc HF equat* ons for the 
system Is very close to or?o of r ho direct product HF states 
of the noninteracting impurities* The energy shift of these 
states due to interaction is calculated tc leading order in A 
(actually tc order A 2 ). For a single impurity, *"he HF states 
can be labelled + , + or 0, i*e., the two degenerate magnetic 

states and the nonmagnetic state. All three are locally stable 
in a certain range of (the mixed, valence regime). For two 

non Interacting impurities therefore, the HF stages are direct 
product states denoted by ( era > , jOO > ■, | a - cr> •, I 0 er> and 

| cO > where the first letter within the kr-t denotes the spin 
state of* impurity l and the second that of impurity 2. Each of 
the above product states, with the exception of the second, is 
doubly degenerate* The first represents magnetic HF states at 
both sites with spins aligned parallel* the second corresponds to 
nonmagnetic HF states at both sites, e^c. In the absence of 
interaction between impurities, the states with parallel and 
anti parallel alignment are also degenerate-. We now calculate the 
energy shift of these states due tc int-oracrion , to order A 2 • 
.This Is In the spirit of the usual perturbation theory, i.e., 

'the effects, due tc the change In HF states caused by A \ and 
the direct effect of A are properly included-; We take up the 
various configurations one by one. 
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Section 3b j TWO ANTIPARALLEL MOMENTS 

In the anti parallel configuration, by symme 
So we define quantities r. and d , 


n 


a 




a ~a 
ry, n 2 * n 2 „ 


• » • (34 ) 


d 



A 



y * r 



»** ( 35 ) 


Equations (29) and (30) th--r reduce tc 


1 / TC 

n - 2^2 " 


tan" 


e d +(U‘2C)n+V 


1 

2 


tan 


-i 


j-i (U-2c)n-7 


4«« 


(36) 


V = 


U_ 

2 % 


( tail 


-1 £ a +(U-2C)n+V 


e +(U-2c)n-V 

, Cl n 

•can — ) 

A 


(37) 


Equation (31) gives, when (35) is substituted, 

V = V( U 2 cl 2 + X 2 ) ... (38) 

Define c i as, 

= £ d + cu - 2C)n + V ii. (39) 


Then, from (36) and (37) we get 

e+ =£ d + dS, ... tail"" 1 -“£) - ~(J -tan” 1 -= ) ... (40) 

Equations (40) are the same as HF self consistent equations for 
a single impurity, equations (2*) and (3*) of section 2(b). We 
must therefore have s + = x + , where (x + , x_) is the magnetic 
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solution to C2‘) and (3 l ); the nonmagnetic solution ( e + = e -) 
is not permissible since then we would have d = 0 or n^ = 
nx° = n| which is contrary to the supposition of antiparallel 
moments. New we compare the energy of the two Interacting, ions 
written in terms of . + , with that of two nonmterectlng iooe 
written In terms of x + . The energy of two noninteracting ions, 
from equation (11.10) of appendix I is, 

F ni = -2 / dZ-'z-n' ( z) -2 / dZ*z*n'(z) + x + n(x_,_)+x n(x_) 

-2U n(x + )n(x)+C {n(x + )+n(xJ } 2 ... (41) 

The energy of two interacting ions in the HFA is given in 
Appendix III and becomes, after substitutions from equations 
(34), (35) and (40) for antiparallel spins, 


a. p _ 2 a 


= 2A[ / 

TC L 


0 


ede 


I) ( e 


^ 2 , 2 

-£ + ) + A 


+ / 


(a- 


ids, _ 

) 2 +a 2 


P ] -2U(n 2 -d 2 ) + 40n 2 


(42) 


The first term in (42), that within square brackets, can be 
reduced tc the same form as the sum of first four terms in (41); 
and since x + = e + , these cancel wh^n F n 1 is substracted 

from F 61 * 5, . So we have the energy of interact! on of two 
antiparallel impurities 

Bint = * Fnl = -2U(n 2 - d 2 ) + 4Cn + 2Un(x 4 .) n(x_) - 

2 

- c n(x+.). + n(x_) ... (43) 

In equation (43) the last two terms can be expressed in terms 
of X+ through the self consistent equations (2 r ) and (3 T ) , and 
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the first two terms in e + through equations (38) and (39)* 
Then most of the terms cancel and we get 


E 


cr-cr 


2 y 

U 


J int tT" ...(44) 

The physical origin of this term is well known and corresponds 

to a coupling via an intermediate sta^e whe^e both + and + 

spins are at site 1 or at 2; the intermediate state thus has an 

extra energy U , and hence the energy denominator U in (44). 

This result is exact within the HFA and differs from Alexander 

and Anderson's result by terms oC ~ 2 ~" ) + higher order in A ; 

£ d 


Section 3c : TWO PARALLEL MOMENTS 


When the two ions are In parallel configuration, one has. 


a a a 

n l = n 2 = n 



v a = X 

• • * (4 5 ) 


The self consistent equations (29) and ( 30 ) then become 
n 0 = 1/2 -(l/2ic) tan -1 [{e d + (U-C)n" 0 - Cn 0 +X }/a] 

-(1/2tc) tan" 1 [{e 6 + (U-C)n" 0 - Cn 0 -X}/a] (46) 

where a = + 

If we put In (46) 


e d + (U - C) n~ r - C n* 


then equations (46) give us, in lowest order In x , 


( 47 ) 


= £ , + 2 =C{" 


tan' 


-1 


'+ 


Ae 




+ 2X‘ 


& mm — C. n « ~ l V\ — UW XX » C. /V r-\ O O 

+ d 71 - 2 A ul + a V 
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l<! - tan “ x -f + 2x2 


A £ — 
+ 


2\ 2 


(ef + O 


... (48) 


The HP energy given by e.g. (AIII'. 2) of appendix III can be 


written in terms of e + as, 


7 parll 


2( (e + ~e d ) (e_-e e )>/(U-2C)-0(e.,.-e ) 2 /{b(U-2C) > 


CL 

& I (e 2 ^ 2 )- 1 + i.5 l !° 

Tl CT 1 % L 1 


do 


2 — t 6 = ± ••• (49) 

o ~'D (e-e) +A 

a • . . 


Comparing this with F" * we find the interaction energy for 
two parallel moments as, 

E int £ FPSr11 - F r ‘ a = -2[6e + (e_-e d )+6e_(e + -e d )]/(U-2C) 


2C(e + -e_) 

UCU-2C1 ( 6e + -6e_) 





-1 


^ I 6 E„ / EdE i(( £ - E J 2 + fi 2 > 

it o ds a 


-1 


a = 


... (50) 


where we have put Ss + = e + - x + , x + being magnetic 

solutions for a nomnteracti ng Impurity in HFA. From (48) we 

2 

can see that e + - x + a. o( X ). Therefore, in order, to 
evaluate (50) to lowest order in X , we replace £ + in (50) 


by x + , and obtain, 

%nt =- 2 l 6£ E ( V £ d )/|n - 2 C) 


20(x + -x_) ( 6e + - 6 e_) 

" ~ ~uTu- 2 cy 
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\ 2 a 

TC 


l uWr 1 


2_ 

n; 


I 


6e a / eds {(e-x^+A 2 }" 1 


0 = + ... (51) 

The sum of last two terms in (51) is equal to, 

2 <S£ + . n(x + ) + 2 5^, n(x_) 

Writing this in terms of x + using nonint^ractirg It HF self 
consistency conditions, and substituting in (51) we find that 
the only syrvivlng terras give 


OG o A , 1 1 o 

Eint = “ X ” ("2 T + "2 ~'2‘ ) - “ X (p + Jr p - } ••• (52) 

where p + are density of d -states at Fermi level when the 
levels are x + * 

The expression (52) has the sam** form as that found by 
Alexander and Andersen and shows that the exchange between two 
parallel spins proceeds necessarily through the d-conductlon 
mixing* 

When both ions are in the nonmagnetic HF states the shift 
is found exactly in the same way to be 

Sint = - 2 * 2 p 0 • • • < S3 > 

p q being the density of states at Fermi level for the single 
ion nonmagnetic HF state* 

We now proceed to find the interaction energy when the two 
ions are in different valence states. 
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Section 3d : TWO IONS IN DIFFERENT VALENCE STATES 

The general expression for HF energy of twc interacting 
ions, found in appendix III is, 


F C2) = - U(n]_+ n x + + n 2 +n 2 t) +£ (n^ + d x + ) + (r^ + + n 2 + ) 2 + 


<0 o; o 


+ - l f- 

% M ( 


£-£ 


a 


-(U~C)n a + Cn" a ~aY°') 2 + A ‘ 


a = 1,2 

Let us first examine V°' . By definition 


... (54) 


V = ^[{(U-C)! 0 - Ccf 0 ’} 2 + * 2 ] 


(55) 


By definition of d , 

(U-C) - C d” CT = (U-C) (ni - n 2 ) - C (n£° - n 2 ^ ) ... (56) 


Suppose ion 1 is In the nonmagnetic low at valence state and 2 
is in one cf the magnetic states. Then (56) gives, if HF levels 
remained unshifted from their nonin;«r acting values, 

(U-Od 0 - C d ^ = (U-C) (n(x 0 ) - n(x + )} - C {n(x c ) - n(x_)} 

where x 0 is the nonmagnetic HF level and x + arc ~he magnetic 
levels m the VF regime- for a single nenin ter acting ion. With 
the help of self consistent equations for e single ion, the 
above quantity can be written as 

(U-2C) n(x Q ) - {(U-C) n(x + ) - C n(x_) > = x 0 - x_ ... (57) 

In the HF analysis of one impurity wa showed that when 
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C» it A , t,hp stable- magnetic solutions are in T go^ f magnetic 
regime, i.e.,- in the VF regime around e d -0 , x_ 1 £ d - C and 
x ^ e d - C + U.' Also,’ in the VF regime, x 0 *e d • So the non- 
magnetic level x 0 does noe come closer than C to either of the 
magnetic levels. So, for sufficiently la^ge G and sufficiently 
small i , we can expand (55) in powers of X , 

X 2 


y T a - ('U-O d 0 - Cd u + 7 a „ cr , 

v - iu ; - 2{(U-C)d -Cd c ) 

Substituting (58) into (54) w c g^t 


-i- 0 (* 4 ) 


... (58) 


f ( 2) = _ u n i + - +n 2 * n 2* + C[ (nt-+nt -) 2 +(n^ 2 j/2 

e de 


+ l 4- ! 


0 


a 9 1 


TC 


D [£-£ d -(U-C)nJ-Cn?+(i -|) X 2 /{(U-C) d 0 ~Cd 0 >] 2 +A 


O n - 0-1 -T 2 , 2- 


i = 1,2 


... (59) 

Tho energy expression for two nenintoreeting ImpuriMss is the 
88Be ea (59), except for the absence of ‘he terms proportional 
to 1 in the denominators of Integrands within the square 
brackets, and that, the n ± 's have ‘heir noninteracting HP values. 

So, we can straightaway write the interaction energy as, 


B int 


_ f C 2) . p 11 * 1 - 


= - U l (&vn?)+C I (I 6 n?)(E n?') 


i o a 


i a 


+ l [{ -(U-C) 5n a + Con" 0 - } 

a t X (U-C)d 0 -Cd 0 


• f a de ^ '{(e-Oi 0 ) 2 + a" } j 

" I) i = 1, 2 


,2, -1-, 


... ( 60 ) 
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To fine? the 1 lowest order correction to F, we substitute in (60) 

0 0 

the noninteracting impurity values for £ i and nj and obtain, 


^ 2 ^ 2 -f 

— Axl_ + 0 ( on • 6n ) + 0(X 2 - ed) 

lnT (U~C) d+-Cd + (U'-C)d + ~ Cd i 

... ( 61 ) 

The last two terms in (61) can be shewn to be weaker than x 2 . 
If we take nj = n£ a . = n 2 (x Q ) and n + = n(x + ), n 2 + = n(x_), 

then 

^ n(x 0 ) - n(x+) ^ _ n( x 0 )_- n(x_) 


Using these and the single ion self consistency conditions, we 
get 


A 2 [n(x 0 )( 7 i- 

'O 


1 


c o" z + 


-) 


n(x ) 


-cr ^ , -\r 
A 0 


i 

X 0“ X -!- ' J 


... ( 62 ) 


In the MV regime around = 0, as we have already shown, 
x + ^ - C + U and n(x + )<u a/( e -C+U). Theref ore we emit 

terms involving x + 4 n (62) and have 


Oct 

E int 



n(x + ) - n(x„) 
x c - x_ 


* • • ( 63 ) 


If x 0 > x_ , n(x 0 ) < n(x_). 

Thus the interaction energy (63) is always negative, and is 

X 2 

roughly 'v, - — 
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Section 3e ;A BRIEF DISCUSSION OF RESULTS OBTAINED FOR A PAIR 


OF INTERACTING IMPURITIES .IN HFA 


The results obtained in the above subsecti erjs are 

<30 2 / N 

Smt = “ x (p+ + P J 

00 p 

E*„+ = ~2X^ p 

1 nt o 


E 


C-0 

int 


E. 


Co 


int 


= ~2X 2 /U 



n(x Q ) - n(:: Q ) 



n 

A 


2 


C 


upto lowest order in x . 

The first two of th-s-o, the interaction energy for 

parallel configurations, ar« dependent on itinerancy (the 

parameter a ) and becomes string ^r as tb- relevant HF level 

moves closer *o "he Fermi surface or p increases. As pointed. 

out by Alexander and Anderson, this is & Zr-ner typo of exchange. 
o-o 

Eiint does not depend on a . Hcwover, in perturbation theory 

a a 

one expects to have further 'x-rms as in E-fp* since all the 
mechanisms for parallel spins ar^ applicable h~r". However, as 
we have mentioned in section 3b, our result is exact within the 
HFA and differs from Alexarr^r arid Ana^rsonfs in nonleading 
terms. For U = « , the interaction energy for anti parallel 

configurations vanishes and "he parallel configuration is 
favoured throughout "■■he range of e ^ » Bu* for finite U, anti- 
parallel configuration is favour^ when 


X 

u" 


X 2 (p + + p 



or, roughly, when 
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(e 


d 


. a / jl , — ~~ < 1/u 
0) d + 


In the parameter range cf our interest, 


or when |e - C | <V 


AU 


In the mixed valence regime cf our err" er^-s- , - G < - C < 

- 2 , as found in section 2. So, if c<V~~ , the 

antiparallel configuration is favoured throughout the MV 
regime. If C > V~~ , the antiparallel configuration is 

favoured when C -V ( au/tu) < s d < - (upper portion of MV 

regime), and parallel configuration when e a < C - ^(AU/n;) 

When £.< c/2, p ip, • Then r, ho state !oo> has 

0. v 7 0 ~r 

energy lower then either parallel or arrl pa rall-1 state. 

The energy shifts obtained above for sp-clfi^d individual 
impurity HF configuration car. be viewed as diagonal terms of 
an interaction Hamiltonian, e.g., E^t = < aa l H i n t I 00 > ‘ 

Th^re is an important off diagonal term left out in the 
above discussion. This is * N h' 5 coupling between the neninter- 
acting HF states | Oa> and [ <jO > . In the perturbation theory, 
the X dependent term in ' he Hamiltonian removes the degene- 
racy between these **wo states and introduces an energy differ- 
ence of 2 X between the symmetric and antisymmetric combina- 
tion states. However, in the HFA this effect is not present 
and we must put <Oa|H iirt |a Q > * x 


os the most Important 
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off-diagonal term, causing transitioiEb ''tween the two HF 
states. In the next chapter we make use of such effective 
Interaction Hamiltonian matrix elements t 3 discuss a lattice 
of such ions. 



CHAPTER III 


THE MIXED VALENCE LATTICE 

We have soon .in Cfcap-^r II the* for a C ^ 0 Infers c-n model 
Impurity, th®re are three locally stable s*.'^s for a certain 
range of values (mixe-- valence regions* actually there are 

two mixed valence regions, one abound i - U and another 
around * 0* W® are in" ~r -st?-' in •-he latter). Two of these, 

corresponding to magnetic solutions ( < hhq > <n^_ c -> % 0 ) 

are degenerate, while th® third low valence nonmagnetic state 
( < n rfcr > = < n^., 0 > <v0) has an energy which may be lcw®r or 
higher than that of the former two depending on e d . Further, 
the energy shifts due to interaction betwe-n ■ K wc such 'mixed 
valence' impurities have been calculate*. Thes® results are most 
simply expressed in a psuedospin (S = 1) model. m this form, 
generalization to a mixed valence lattice is imm -fJ.fte , and one 
can discuss the thermodynamic properties of a mixed valence 
lattice, especially oho- question of possible phases of the 
system. This is done below. 

Section 1 i THE MODEL 
Section le : SINGLE IMPURITY 

Consider first a single impurity at site, i . The ~hree 
locally stable HF states are denoted |+> , ]+> 9 |0> a nd 

have energies E 0 , E 0 and E 0 


respectively. E a and E 0 



can be calculated i n terms of the model parameters A , U , 

C (see Chapter II, sec. 2). This manifold .Is i ! Fp»' , --'S?.n i ‘i?d by the 
states of an S = 1 system, with a single Jon anisotropy term, 


3 .e. , 

H = ( E a 
The states S z 
magnetic HP states. 


i 2 

- E 0 ) S z + E 0 ... (1) 

= + 1 are d^gen^rate and represent the 
S z =0 represents the nonmagnetic state. 


A Hamiltonian of the form (1) with acV'i^ tens 1 off diagonal 

1 14 

terms of the type * S z has been used by de chafcel as a. model 
for a mixed valence Impurity. Such a term leads *-o the magnetic 
susceptibility being always finite at T = 0. However, if is well 
known that the effect of terms beyond the HFA is to cause tran- 
sitions between the HP eigenstates wi 4 h long dime transients 
(Anderson and Yuval 9 , Hamano 8 j'this thesis, Chapter iv) . As a 
consequence th^r° is Hondo quenching of ■'‘he moment, at low tem- 
peratures. Their effect, carro 4- be simulated by a simple off 
diagonal term such as the above. 


Section lb : TWO IMPURITIES 

Consider the following model Hamiltonian for two inter- 
acting S = 1 systems at si^c-s i and j , 

■f i i ^ n’2 4 2 \2 

H = J 13 S z S^ + Bjj S z + ii(Sj + ) + M ... (2) 

This Hamiltonian is diagonal in the S* , S*j space, and the 
eigenstates are in one tc one correspondence with the two 
impurity HF eigenstates discussed in II. 3. Further, we see 


& w- ! , , 
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that 

the interaction energy shifts 

obtained 

th~re are. reproduced 

exactly by (2) , if 




= * - 421 ) 


♦ .* (3a) 

B i0 

= £ ( B i n * + E int ) - 2B i°rSt, * 

w 

Hot 

i.. (3b) 

A 

aU 00 

= % - s c + 3 int “ %nt. 


**# (3c) 

M 

00 

- 2E C + Bj nt 


(3d) 


The interaction x , being th- 

> H -f v&ct 

overlap b ?+, we s n f 


orbitals on d iff pr^n f ; plus hopping via conduction elec- 

trons, is very small. We therefore consider only rr rest 
.neighbour interact,! one. 

iis wo have pointer? cut, in Chapter II, Sec, 3, the inter- 
action X has off diagonal m? ^ -rf x el--m'r;*"S In the HF basis, 
between scares | Ou > a :o [o'O’ « This is easily included in the 
pseudospin model through " h ^^m 

= X Sj S* ... (4) 

Sine'- the ax 4 s of spi n GU<snti za-b 4 on is 2*rb rfc^ary and the 
system is actually rotati onally Invariant ("hough no": 0 HFA) , 
th^re ought to, be a " ransv ■ "•? spin coupling in s-**]** to 
the long! : 'udi nal or isfng coupling bo-w^^o. *h- spine. This 
coupling flips *-h<- HF state at site i r fom |+> -o j+> (say) 
and simultaneously the HF state at site o' from | +> r Q |t> • 
Its magnitude is easily obtained from *he condition of 
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rotational invariance. Combining all th^se terms, the results 
obtained in the HFa for two interacting MV impurities can be 
described by the model Hamiltonian, 

4 ' .2 .2 .2 .2 * , 

H = J 1;j S z Sg + B i( j S z S z + ACS* + $z ) + Xs x S x + M 

. ••• ( 5 ) 

where the coupling parameters &m expressible in ^crms of the 
basic Andersen mod^l parameters , & , U , G , X * 

An objection can be made that by introducing (4) we 
completely include the eff • ct of X , and hence cth^r terms 
are superfluous. This is not' so, Th<~ &n + 4ferrcmegnrii c coupling 
arises from Intermediate s* etes of type = nj + = 1, 

not included in the thr^e l^val manifold b-^ing used hero. The 
ferromagnetic and other couplings (viz,, SB » <5E 0C f » 
6B 0 o) are results of itinerancy or mixing with conduction 
electron states. That also is a feature net included in the 
three- level scheme. 

The Hamiltonian (5) does not have two important features 
and because of this it is a poor model f or explicit dj.scussl.on 
of valence transitions in *he concentrated (lattice) system-. 
These missing features are the following : 

(i) In mixed valence systems the transition r c a low 
valence state takes place with promotion of the extra electron 
to the conduction band-. For a mlx^d valence impurity this 
leads to nothing new-. However, in a concentrated system, th^ 
fermi energy becomes depend ^nt on ^ho valence fraction. In our 
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model, where the fermi energy is always kept zero, such an 

electron count effect‘d c& i be phenomenologically taken into 

5 2 

account by making e a function cf < S_ > or the valence 

U. Zi 

fraction. It is generally believed to have • K he effect of 
pinning the d level + o "“he fermi level. The functional 
form is to be computed from knowledge of density cf conduction 
band states etc. 

(ii) as pointed out be Chui and Andersen-^ ions cf differ- 
ent radii in an elastic medium (the solid) cause long range 
strain effects due tc size mismatch. Presumably, this can be 
partially taken into account by an additional long range 
coupling term S * 2 S ^ 2 . 

Because of these shortcomings, we do not use the above 
Hamiltonian to explicitly discuss valence transitions as are 
observed in S m chalcogenid es , Ce , Ce x Th2_ x etc. However, 
the model is quite suitable for a ciscussion of the phases 
possible in a mixe^ valence lattice. We not®, though, that an 
important phenomenon, namely the Kcooc off ?c A , is missing in 
this HF^ modelling cf the system. W;- discuss the consequences 
of this omission later, end in Chapter V suggest a scheme which 

t 

includes the Kendo effect as w<=ll. 

A Hamiltonian of the form 
2 

H = 6Sj + m Sj + J Sj . S j 
1R 

has been used by Ghatak to discuss the valence transition. 

X 

This has a single ion transverse spin term «Sj whose micros- 
copic origin is unclear. It r’ms not have the longitudinal and 
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transverse near neighbour sp-*n - spin couplings which wp have 
found In Chapter II. Fu' TV "hrr > -h-' twc f«a**ur°s mentioned above 
are also massing. It c a -net ^h'-r^fere be a gr-v* model for dis- 
cussion of valence •-ransi-^ or for various possible phases. 


Section 2 s POSSIBLE PHASES OF THE SY "’"FM 
Corsid er th^ Ha.nl l" cni an , 

* 2 1 i v f2 «2 

H = * I S£ + M + Z J» j S z Sg + Z Bj j S 2 Si + 

1 i,3 _ 1,3 

+ x I S* S^. ... (6) 

13 

This represents an S = 1 lattice system with 
(i) single ion anisotropy 
(11) anisotropic exchange 
•(ill) a kind of quadrupolar interaction. 

Such S = 1 mod ils (generally with more symmetric ccupllrgs) 
have been extensively discussed m the 11 *• 


The possible phases depending ers system parameters., are 

i 

(1) Ferromagnetic, viz-., < S 2 > ^ 0 

(il) Antiferromsgn-tic, viz., < S z 3u1:ila '‘ k '' 2CP > 

(ill) and (Iv) Similar phases, with alignments In the xy 

i 

plan-"', e.g-.', < S x > /o 


(v) 


Singlet, i-.e-., a 


< 


a 


sublattice 


i 

phase with < S a > - 0-, 
= 0, a = x, y, z- 
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and < S z > (we have assumed h=re that j is negative, i.e., 

the coupling is ferromagnetic). The Curie temperature is 
obtained by requiring that self consis'-eir 1 ' scluti oos wi^h 
< S% > o exist. In this way one obtains the self consistent 


equations as, 


-2e 
l+2e" 


,(A+E)<sf>) 


z 


cosh<p J < S >) 


< S z > = - — - . : — — sinh( p J <S > ) ... (8) 

z cosh (pj<S>) 
z 

~p(A-r 

o™ — eosh(pJ<S >) C9) 

-P(A+B<3>) 

cosh(GJ<S >) 

Z 


1 + 2e 


where B = I and J = I J* , and wp have emitted 

j (run.) J j(n.n. ) ' J 

the site index. 


£ simple special case arises when B-j j = 0. In that case, 
one finds that the self consistency condition for < S z > , 
e-quation (8) gives the condition for < ^ o to be just 

possible as, 


~2j3J e"P A 
1 + 2e“^ 


= 1 


CIO) 


This equation determines T c , which can also be see: o from the 
usual susceptibility ?xp~esei on in MF* - 


X” 1 (T) = X” 1 (T) + J 

-1 


(ID 


where x Q (T) is the susceptibility of a single ncnint^racting 

s ca i 

2J3£ 


ion and is in this case 

■pA 

M 

X 0 


1+2 6 


-pX 


( 12 ) 
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So condition for the susceptibility to blow up is the 
same as eq. CIO) above. From this equa tl.cn , it Is clec.r that 


for 1 < 0 

k B T c 1 ! J 1 *“ 

the well known Cuvie-W- : iss result-. For L > 0, the transition 


temperature decreases exp on -nt islly , satisfying the relation, 

-a./(k s x c ) . ( 14 ) 


JD 


hi e 


It is well known , however, Ohs'*" for a. > k^. the transition 
is first order, the point A. c -, T c C-A^) being the well ±cnown 
tricr.it! cal point Further-, for > -.463 |j| the ground 

state is nonmagnetic (a singlet or < S Z > =0 state) , since in 
that case no real solution exists to the equation (10) which 
determines Tc. 


In our calculations-, th-~ results of which are shewn in 

3 

fig-. 8, we find that generally £ * A a tid •] Jl^lO . 

(for the chosen parameters II, G and X ) unless we am very 
close to the HF phase transition line where 1^0. Thus-, the 
results which are crude solutions o equations (13) and (14) and 
basically show that as e d increases within the MV regime, 
starting from its lower r-p&ch, o q. slowly increeses j o a maximum 
around the HF phase transit"! err, and “h-n falls off exponentially 
over a range 10” 3 a Th® assumption of B = 0 is not un- 

realistic since we find B * 1Q~ 3 &«A ever most of the range'. 

In any case, it can only shift the peak in T c * ^nd will not 
change its character in any significant way-. 
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Section 2b : THE VALENCE FLUCTUATION PHASE, < S x > / V 

Tc investigate the < S x > ordering In MFT , we consider the 
mean f i e Id Hami It oni an , 

4 2 4 

H = J J.±‘ Sg. + • * « (15) 

i i 

where, A* = A + B < S z 2 > 
and \ x - I *i j < s x > • 

j 


Thus, we are looking for a kind cf ordering wh-re < S z > = 0, 
< S x > ^ 0 . Within the mean Held theory re have the usual 
expression for susceptibility, 


-1 


-1 


X x 


+ 


(16) 


where x° is the single(nonintorscting) ion x susceptibility 
x 

and Is found tc be 


. o 
'x 


4 -l 


4_ 1 

A* 


l-i- 2e 


■PA’ 

=PT 


k_ 

8 


[1 3 coth £——] 


So the solution + c the equation 


-1 


= 0 is 


... (17) 


(18) 


k T 
B c 


-1 


2 tan /A' - 1) 

( 2 


... (19) 


From the above expression it is clear that when 2A <A'(A*>0) 
the transition temperature is nonexistent or negative, so that 
the system oo^s not order in this particular mode. 
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Fig* 8 shews the results of calculation cf T c from (19). 

It is nearly coaster.* over 'he MV regime and is ^ % A as 

d 

expected from (19). 

Section 3 : COMPETITION BETWEEN PHASES 

We new discuss which of -he possible phases ~he system 
will actually have a* a given temperature end for a given set of 
system parameters. The problem can be approach"-'’ f rr. m the high 
temperature side (com pa riser, cf T c *s) an? from the low tempera- 
ture side (ground state). 

The comparison cf T c ’s for system partmer^-s In the MV 
range are shown In fig. 8. For simplicity, w e consider C < ]fU A 
so that only Ferromagnetic ordering is possible. We find that 
ov-r almost the whole range, *he <S X > ordering dominates. The 
possibility of <S Z > ordering arises only at very low tempera- 
tures. We show b a lew that such temperatures are too close +o the 
Kendo temperature to be meaningful* since* as ws have already 
mentioned * the theory does not held an *hese Jr "mp<= natures and 
below. 

In fig. 8, we show <-hi varie*l on of critical temperatures 
f or < S x > and < S z > orderings * The parameters ch 'sen are 
U = 200 A , C = 8 A , A = ,25 A , Computer calculations show 
that the MV regime extends from ^ ,7 tc * 6. When «w 6, 

the magnetic HF solution is (x, y) «v (- 1.5 A , 160 A ) which 
corresponds to a Kendo temperature, Tg; ^ * 


15 a .At %he 
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other end of the MV" regime, th<= spin fluctuation lifetime is 

expected tc be v • These tempers tu res are below the x- 

x 

ordering temperatures (Tc ) but wp-11 above the z~ ordering 
temperatures. Thus, if at all any ordering occurs, it will be 
the valence fluctuation type of phase. If X is too small, say 
by another order of magnitude than that considered here, this 
phase is also unlikely tc occur. 



CHAPTER IV 


BEYOND THE HARTREE EOCS APPROXIMA II ON 
THE PATH INTEGRAL METHOD POP THE A-TDERb OIT MODEL 


Section 1 : INTRODUCTION 


The HP A discussed in the previous two chapters suffers 
from a serious shortcoming. Per a single impurity ..described 
by the Anderson model (with U » tca ) , it predicts a magnetic 
ground state (for < o) , whereas experimental evidence avail* 

» able in abundance shows that the state changes over smoothly ’ 

' into a nonmagnetic ground state well below the Hondo tempera- 
ture Tj£ • Thus the conclusions derived in the previous chapter 

X 

are unreliable if, e.g., the transition temperatures T c or T c 
are of order Tr • In this and the succeeding chapter we attempt 
to include the Kendo phenomenon in cur theoretical description 
cf dilute systems described by the anderson model. We largely 
concentrate on the case C = 0. 


\ 

Schricffer and Wolff ^ (SW) showed in perturbation theory 


that the magnetic properties of the Anderson model with < 0 


and | £^| >>■ A , are the same as those of the Kondc model, 

I \lr» C V^r C lr! « *** (1) 


H = ko ° ko ° k ‘ I " J ^ + kDcr '^'“ka“k'a 


With j !2|v kd l 2 (i- _ i 


'd 


e d + U 


) 


m m m 


( 2 ) 
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v kk’ ^ v ^' 


ft 


kd 


+ 


_ 1 _ 

* 4 - 


fr) 


V.'. (3) 


Here, S is the spin of the local moment and s is the 
conduction electron spin density et the local moment site'. 


It Is w^ii known that a Kendo impurity described by eq.(l) 

has a nonmagnetic ground state. Apprcr-mately, x + Tg) 

-1/Jp 

where T = D e ff '{ J p e for J p « 1, D eff being an effec- 

tive bandwidth* The correspondence between And-rsT, and Kondo 
models is strengthened by the results of Krishnamurthy, Wilkins 
and Wilson^ (KWW) . Doing a renormalization group calculation 
which amounts to a numerical nonp^rturbative calculation of 
energy levels of the system, they shewed that the symmetric 
And erson model (e d = - 2 . and in (1) and (3), V - 0) is 
identical with the Kendo model in Its thermal magnetic 
properties, with J given by the SW value, with D e ff = U/12. 

For systems of interest to us, viz., U » rc A , -*■» 0 t 

the SW transf ormation is of limited value, since Jeff *► ® 
as 0; this explosive situation arising because of the 

breakdown of the perturbative basis of SW transforms ti on 
as o. The question then is whether the connection between 

Anderson model and Kondc model (-for magnetic properties of the 
former system) breaks down or is modified in a definite way^ 

KWW have studied this question in their numerical renorms liza«* 
ticn group approach and have found deviations from Kondo like 
behaviour for |a d j < A .We study the same problem here in 


i 
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an analytical, albeit more crudely approximate, framework* 

The method we use is a generalization that due to 
Hamann® for the symmetric Anderson model* This method has as its 
special feature, the starting point of Hartree-Fcck states 
rather than the ground states of an unperturbed mod el which are 
then perturb at ively expanded either in terms of A |e^| or U/ A . 
Such perturbation expansions would be useless for the range of 
parameters of interest to us. Hamann's method then proceeds to 
describe the system in terms of T paths* or fluctuations between 
the HF states arising out of processes beyond the HFA. in his 
treatment he has chosen only special types of paths - linear 
hops between the HF states - which correspond to sudden (though 
not discrete) changes in the effective d electron HF potential. 
Nczieres and De-Dominicis®® have shown that such changes excite 
low energy electron-hole pairs in the conduction electron system 
which are responsible for Kendo divergences. With ^his approxi- 
mation, Hamann has shown that the partition function for the 
symmetric Anderson model is the same as that cf the Kendo model 
as obtained by Anderson and Yuval 9 . This approach shews that- 
because cf the near orthogonality of the two HF states, the 
” communication* 1 between them via excited electron hole states 
is small and Tg is, consequently low. Hamann obtained the Kendo 
model parameters for a given set of values of U and A (he 
has taken e d = - U / 2) . The parameters are, 

J p =• 4A/U ... (4) 

D e ff — U / 12 ... (5) 
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agreeing completely with the later results of KWW. We extend 
here this approach tc the somewhat mere complicated asymmetric 
case • 


Section 2 : THE PATH INTEGRAL FORMULATION 


In this section we closely follow Hamann and obtain a 
formal expression for the partition function, rhe evaluation of 
which is done in the next section. 

The Anderson model is broken up as 
H = H 0 + E 1 

where H, = J ^ n k(J + J e da n aa + Un^n ... (6) 

ICO o 

anfl % = V 2 ( c ka c da + °te c ka> ' — < 7 > 

The partition function can be written as, 

Z = Z Q < Tj exp ( ( t) ) do > ... (8) 

where the quantity after Z Q is an abbreviated notation for 
the usual perturbation expansion, T t denoting time ordering, 
and 

H % -Hr 

% ( t ) = e 0 H-l e 0 ... C9) 


when right hand, side of (8) is expander,, we have sums of time 
ordered products of the form ... )%( to) 

integrated over the t's * The operators H^ change the d 
electron occupation numbers n ^ , but H 0 conserves them* 

therefore, within the trace (8), exponential factors arising 
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from (9) contribute numerical factors cf the f crm e“ Url d + n &+ x * 
from H 0 , with n^ now denoting the eigenvalue of the corres- 
ponding operator. The product + n^ + can be written as 

n<3 1 034 . = -£• [ (nfl + + n^ + - (ndt ** A3+ ] • • • CIO) 

then the above exponential factors Involving ryt n^+ can be 
written out using Gauss integral, formula, 

e 8,2 = f dx e “^ x 2 “ 2 f^ ax 


After carrying cut these steps, we get the expression for Z as, 

, , -/ p tt + t )clT 

Z = Z /5X/6Y Z + (X,Y)Z + (X,Y)e 0 ■ U U (11) 


In (11) , Z + and Z + arc- partition functions for + and + 
spin d electrons with d energies, e d + X + 1Y, with the 
+(-) sign for + C + ) electron, mixing with conduction electrons 
via the mixing given by H]_ (equation (7)). The U term has been 
removed by the above mentioned steps and results in the general 


time dependent fields X and Y . Thus X and Y are general 
functions of v . 


To evaluate Z (expression (11)) for a general path, we 
look first Into Z (X, Y) , which is given by 

Z a (X, Y) = e Ca ... (12) 

where C a is the well known closed loop sum. This can easily 
be shown to be equal to (e.g., see Hamann) 



• * 9 


(13) 
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where Av = A(iT(T) + 'aX(x)) is the external potential 

(A) 

acting on a d electron of spin a • if the Green’s 

Turction for the cUelectr on in this potential. 


Formally, the Green's function for an arbitrary time 
dependent potential can be fourr from uyson’s equation, 

G = G° + G° VG 


(14) 


■where G° is the Green’s function in th" absence of ary time 
dependent potential, i.e., when the d electron is In the level 
% (say) and mixes with the conduction band as usual tn the one 
body hf theory. G^ 0 ( r ) is the Fourier transform of G^. ( a ) = 
(e - E^q- - i A ) 1 , and is found tc be, in the long time 


approximation, 


o 

G d 0 ( T ) - 




... (15) 


o 

The Dyson equation (14) with G^ 0 given by (15) and with Ed = 0 
has been solved by Hamaen, using Muskhellvilis 2 ^ method for 


singular integral equatl ~ns . Wo also put E^j = 0 so that the 
instantaneous potentials for a d electron of spin a is, 


V 0 = £ d + iT(x) + aX{%) ... (16) 

Hamsnn’s result for G^ (t) is (we have dropped The spin indices 
which are to be substituted back wh<=n necessary), 


G d (x > T ’ ) - ®a(T,T', 5 (x)) + 4- It - £ ~ ra - [-4— 


1 


V 

Ji. 


A' ( t) A*™( t) 

[ A ( t ) - a|(t ! ) - A ~( x j ) j 


] - 


, ”1“ , 


T — T 


(17) 



whero 


S( x ) is a dimensl onless quantity defined as 


C^(t) = V^(x)/A » * * (18) 

and A"(x) = exp (+ir|( x) a(x) ••• (19) 

where A ( x ) = exp [~ f Sill. llll-Ll ] . . . (20) 

L % ' x ' -x J 

with rj ( x ) = - tan“^ £(x) • •• (21) 

After having four d G dc >( T » o'), the next step towards 
evaluating Z Q is to fin^ ""lie linked cluster sum d 3 given 
by (13). Hamann has shown that this can be broken up into an 
adiabatic and a transient part corresponding to adiabatic and 
transient parts of G$j — the first and second terms of r.h.s. 
of (17). respectively. The adiabatic part of C a , which we 
denote by.’^ 

M a = - i f P 5 a (t)[f - tan- 1 ( m (l+e|(x>)3 - 

0 <* ... ( 22 ) 

The transient part of C a is denoted by B a and is obtained 
by substituting into (13), 



1 

% A£( x) 


[Sin r](x) 


drj_(ji 

dx 


sin T) ( x) 


cos r)(x) 


a 

ax 


In A(x) 


... ( 23 ) 


The first term of this contributes to 


G a s n Integral - 
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This integral vanishes if rj ( p) = n C 0) which must be the 
case for periodicity of Q d . Sc only the second torn of (23) 
gives & nonvanishing cro^r^hut-i on i o B a * ¥e find that 


B = - JL f a ^ ' f AA s.in_2ji(_T)_ . ••• (24) 

2 v T „ T > ■' x 2 dr : 

TC 

The condition tj ( p ) = rj (0) = 0 has boon used m (24). 
Substituting r) (t ) = tan" 1 xs( t) in (24) th* coupling 
constant in-'gral (over X ) becomes, 

/hi misl 

o {1 + X 2 £ 2 (t)> {1+X^C 2 (t [ ) > dT * 

- i n . iCl’l ... (25) 

2 { C 2 ('c)^ 2 (t 1 )> l+? 2 (n') dT ’ 

and putting this in (24) we got 


B = 


_P_ 

2 71* 


dr dr* 
J A_fr * 
0 


r— r 


d€ ( t 1 
dr 5 


„ _i( jj. 


in 

( t' ) 


ehihl 

l + h‘(x' 


(26) 


So, finally, for Z we get, 

Z = Z 0 j 6X/6Y exp { I (.4 - 1 - B ) - A 2 / P dr I £ 2 /U >•«* (27) 

C O c 

where and B are given by the expressions (22) and (26) 

cr 

with appropriate spin indices put in. 


Section 3 : THE CHOICE OF P.-THS AND EVALUATION OF Z 
Section 3a * THE CHOICE OF PATHS 

Before evaluating Z for general time- dependent X and Y 
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wa lock in* 3 it for timr- independent, X end Y , i#e«^ the d 
electrons of opposite spins moving In fields constant in time,' 
which is equivalent to HFA. It, is, how^v^r, not clear rightsway 
which values of X and Y give us exactly the HP solutions • In 
*he spirit of the path integral method, we first write down the 
fre-e energy cf the system a rV find the values of X end Y fer 
which this has minima, with the Intort ion cf then taking into 
account fluctuations around or among these ratline. 

In appendix IV, th-= fr''- 5 energy for cons tent fields is 
written and ex^remizod w.r.t. X and Y. The extremum condi- 
tions arc found to be the seme as s«lf c ns i stent HP equations 
( 2,3 ) of chapter II* Further we show that the conditions for 
the extrema to be minima are the seme as stability conditions 
for HP solutions (equation ( G ) and (7 ) of chapter II). The 
solutions of these HP equations together with the above condi- 
tions, were shown in chapter II to give magnetic solution for 
- U < e ( 3 _ < 0 (U > > ix A), and nonmagnetic enr-s for rest cf the 
range of a^ . Let us denote by X c and Y c the fields which 
minimize the f roe energy. 


Clearly, small fluctuations around the HF minima do not 
significantly affect the ground state which continues to be 
magnetic. However \ transitions between + and + spin HF 
soluti ons, i.e. , between the. broken symmetry solutions, can 
lead to a nonmagnetic or symmetry restored ground state. 
Therefore, like H&mann,. we make the choice cf hopping paths - 
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the 6 electron initially being in -:n<- r f th~ HF stages and a 
sequence of time dependent potential flipping it tc the ether 
HF state, then hack again, and sc on. Such a T-tmr dependent 
potential is shown in fig* 9. Most of the tim'S th« potential has 
the value £-+• either of the HF solutions and '•he time taken tc 
cross over from on- n * the c'her, the hopping time denoted t 
js short compared to time spo"*~ at these. Such paths correspond 
to long time spin fluctue"*! 'ns. For the time being we k^ep the 
shape of the pot -trial within <=>ech hep arbitrary. This, and the 
hopping 5im^, which we as sum- Ic be the sane for «£ ch hep, will 
be determined in a let^r section. 


Section 3b : EVALUATION OF THE, INTEGRALS _* jJJD . 

We first evaluate the JU erm given by eq. (26) - 
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dTdT'. iw li'i 

t--t ! dr* 



5 .( t ) 

_ o 

(77 "-V^'(V7 


In 


1+5 2 a (x) 
1+5 x 1 ) 


Ir. tht' above we first integrate over r by parts and put the 
condition ( p) = 5^(0) = 0. Then , 
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cr 


2ti^ o 
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aten; 


_ a_ 

CLT ' 


In t-t 1 ^r--{ 

r i 't L 




du 1 ? / 


In 


1 +5^(o) 




1 +5^ (r ! ) 


For a hopping path s 0 - ,lh; constant except during the hops. 


“ :k t 

unless t is inside a hd#,* d| 0 (T),/d'C / = 0. Therefore, non- 

van 1 shir g contributions to (28) com." orly when t end t* 


(28) 

So 


both 
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are msi^e hops. If we denote by the ccttributions when 

is within the hop ant o' within the (i ^ j) t ant 

by S 3 the contribution when v ant -a’ are both within the 
i^ 1 hop, then % can be written as-, 


B r 


= i if 
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(29) 


With 

I a . : 
iJ 
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t . +T 
..1 O 
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2 n; 


ut 
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V T o d£ (t) 


ax 


t . 
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da; 1 


■In 


T— T 1 


«<,(■*) 


and , , 

c . -i-T 

1 1 0 

S. = — 5- f dr / dr 


1+5 ( t) 

111 — 5 — } 

1 +^(t’) 


ttT 'CgtTj-S^T') 


2 , 

O' 


a t 


C l ’ r <*5 (t») d 5 (t) 

1 --- — -mlr-T a 


b . 
1 


ax 


In } 

1 + 5 ^(t') 


... ( 30 ) 


C3D 


Th-NSe two terms are known as th^- long and short range terms 
respectively.- In evaluating L 3 j to lowest order, we do not 
need to know the hopping shape or time* In the lowest order we 
simply write (30) as,. 


a 


L i i= ? 


1 , , t ±+x ° f v\> a? (t«) 

In |t . -d . | / dt f J do; 


1 3 

1 


dv 


1 + 5 ^(t) t 

dn ^ } + 0 (•.“;• -r — r) 

1 + 5 J(t') ' b i “ t j' 


d **il. 

dT 5 a (T) " e c (T,) 


.** (32) 


In (32),- the t integral can immediately be -''one and we obtain,, 


, L 'i n ' T o d£ (t ! ) f 
l2 = (-O 3 In It.-tJ./ dr* — 

13 2 % d J - 1 - 


t . 

1 + 5 ‘ 


a.v - [ E 5 : 5 2 (,Ty 


... (33) 


1 -vo 
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where the limns ^ end ? 2 are th- HF values en* P^tain to 5 
within the bracket. The overall factor (-O ' 5 takes cere of the 

, ^ . 4 c •* C In eV ! -0 ltl r ~ c. TO hops Will 

fact that the initial = tc,te is ^2 v -- 1 

be and l-ft 2 . Next we denote K a (^ ’ ) ^7 *' £nd change 

the integration variable to V an' - wn^e (33) aS 


r 2 5 - . 1 . + A 2 n C 2 

2 ln t i“ t o/ ° 4 ^-V 2 1X1 1 + 5~ 2 J C=?. 


i-w" 1 *'" vV«' I 


... (34 


.. -+■ a -m -r. give" In equation (25). 

Using the coupling cor-st^r. .nt^gr^- *- J v » 


the above cs • b^ written as, 


1 5 


_ . in |t . -t . | / &■ S <35 

2 lt J n 5 1 


o' ' " ' i7x 2 'p i+x fc Cn 


ill 1 1 . -t . I / dx { tan' X 5 2 

i j 0 


i 2 *1 > 

: an X 5 3 H i + x 2 5 2 1+ x 2 ? 2 

kl -J- 


/ \ it- j t *1 > — 1 >- \ 2 

_ — in! t^-t^| (tan ? 2 ~ cail 

2 a 


... (35) 


The method, use* by Hamarm in going frun eg-. (34) tc (35) for 
the symmetric case cannot be - nhe asymmetric case, 

whereas "‘he above *rick mak-s 'hi' 5 possible-. 


Now wo take up 


th-' 5 short r£.rge r'^n, eq. (31) 


V T 0 .V T 0 ac (tl 
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Integrating by ports ever % we get 


s a I^I 1 f'O dE - lT,) 
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u i +T o d5 0 (T') 


/ 
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ClT 


dT 1 


l+CU) 

n O n 

In o"* J 


?_ r T 2 

To evaluate (36) we n-^e' tc know g a ( x ) explicitly. So next we 


1+5?aV 


determine the hopping sht pe or th** shape cf the potential 
within each hop.. 


Section 3c : DETEHMI3\UT I ON 0? HOPPING SH^PE , i.ND OP x Q 

Two types of hop are shown in figure 10(a) and (b)#.We 
choose a single hop (a.) rather than an oscillating one (b), the 
reason being as follows.’. 

The short range contribution from a single straight hop,, 
as shewn by Hamsr.n is a factor 0(J x Q ) An oscillating 

path cf the type shown rfr fig#. 10(b) will contribute e 

, jr. 

factor 0 ( ( J x Q ) . The ling range cortributi on ^cos- not 
depend, on x Q or the shape of hops an* will therefore be the 
same for both' type cf paths#. Therefore contribution* t o Z from 
the second type of path (b).ls smaller Then "'het from the type 
(a) by factor o((J x )^), since j x is generally < 1 • and 
therefore can be neglected. So w= do not consider paths with 
many short time oscillations# , v • 

Hamann has chosen a steady straight path by keeping the 
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valence field, Y , constant. For an asymmetric .Anderson model, 
however, the path should be so chosen that as X moves from 
one HF value to the other, Y changes in such a way that the 
free energy is a minimum w.r.t. it (the Y field). In other 
words, we must minimize the free energy for the static path 
w.r.t. Y for a given X. This energy is w-i-ten as (see 
appendix IV) , 


aae 


X 2 Y 2 A r f 

V - u - + u + % ^ f ", o p * 

a (£“£ d - iY + cX) 2 + A 2 


a = 


.. (37) 


Differentiating (37) we get 


9 V 

aY 


2 I 
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+ 


1 r 

— I m 
7 t L 


1 CQXL 


£ d +iY+X 


tan' 


£ d +iY ‘ 


•A 


.. (38) 


When the r.h.s. of (38) is equated to zero, the solution of the 
equation for Y will give an extremum. To check whether this 
extremum is a minimum wo different. late (38) again and obtain 


3 2 V 2 , A/jL , A./tc_ 

3 Y 2 U (e d +iY-hX) 2 +A 2 (e d +iY-X) 2 + a 2 

This quantity is >0 for all real Y. Therefore, all 


paths found by "he equation = 0 

find the paths w? must, solve 


are minima in Y* To 


= “ (it 

n 


tan' 


-1 E a +iT+x 


Tan 


-1 


e / n+iY~X 


(39) 


2iY 
" if" 
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* • • 
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ica) SoMion to eguaiion for minimum enercjy path . 
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A graphical solution to (89) is shewn m fig* ll. The approximate 
solution is, 

Y '^ '4 + I ’i x K 

v ifxl » Ixl <£ d 

W» take X to be varying ltn^'rly. wi~h tin 1 -, vl”Wn a hop, 

i • s » , 

X(t ) = Xg + ( T „ t> ) ... (42) 

T o 1 

within "'■he i ^ hep, if it Is from Xp to X^ ; cth-twtsc. ? n 
(42) X^ an 1 ’ Xg should be In* srch&ngXi • 

Thus, equations (40), (41) and (42) determine *:he hopping 
shape. 

Wh^n this hopping shop- J? substituted into f h^ short 

range term given by eq. (36), many integrals result, which 

have b=an shewn in appendix 7 * It has also been shewn there 

tan" 1 S 2 -taxi" 1 5- 1 2 

that the only r dependent term is 4 {- ~ — — >1 ii|t I. 

To determine r 0 in terms of the model pa ram<-? ers , we 
collect all th-= r - dependent terms in *h«- exponent- of Z 
for a given path’. These contributions a”*', (i) from ‘'he short 
range term mention^ < abov®, an* (11) from ; h<~ In^-'-’g^a *1 on of 
adiabatic plus Gaussian factor terms. Tlrm this t q - 

dependent contribution is minimized w.r.t. t 0 . This is also 
done Jr appendix V and we get 
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... (4 0 ) 
... (41) 



x ^ 
o — 
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l £ a 


(i 
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u 


) 


the ^sult beirg valid for js^l > 


(43) has -'-.ho KWW or Haman'" vt lue for 
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... ( 43 ) 


The result 
U 

" 2 


The rpmeining short range terras ce’frlbut^ to the ampli- 
tude associate-* with etch flip which can be identified with 
J j_ T o • For th- p yranv a r i c case £r A for U » ft A , Haman n has 

estimated the leading term in this integral, &nS has shewn 
■ A • T 

lj_ ^ ij ^ J 2 T o • However, the coefficient is uncertain, 
this being a deficiency of this method which is poor for short 
time properties. We find our integrals +-o be much more compli- 
cated, but for the U «> , lerge negative case. 


J 


X 



Unfortunately no precise statement car be made of the appro, 
ximate size of J X or t 0 for small except tc say 

that the integrals are all firu^e, sc -‘■h-'So quantities stay 
finite as •> o* 


Se chi on 3d : THE PARTITION FUNCTION 

The partition func r I on is the sum cf contributions from 
all possible paths, each pa^h containing 2o heps (n is an in- 
teger). The contribution from a path with heps at , tg..., 
..., is obtained by collecting co^tributi'-ns from the long 
range, short range and adiabatic terms : 



82 


where Z =(J x ) 2n /dt 1 .=d 


n 


exp{ I ( -) 1 “' 5 ( 2-e)ln 

i j 
-1 


tan ^ S 0 - tan 2 

S . £ = 2 ( 2- - 1 - ) 




0 


This form cf th~ pc-r^i-ica func r i on 1? 1 ■- 1 pr* j ca 1 with that 

9 

cf ^he Kendo model obtained by Yuvcl and And •■- rson . The Kendo 
model parametars, 1; , e , etc., car. be obtained as functions 

of Anderson model parameter's , , A , U. 


It is clear that es a d -*■ 0 w« gc into a very strong 

coupling Kendo model. However, effective J z remains finite 
for e d -*■ 0 . This is an percent qualitative conclusion. We 
also assume that isotropy continues sc that, j z . 

However, when J . t is ls-’-ga, flips are frequent, an' 1 the 
above picture based on well separate! hops breaks down. So the 
Anders on- Kendo analogy is, unfortunately, not very accurate 
as e d G . 


Using the above Anderson-Kondc connection, and known 

results for Kondo model, we can calculate T K , x(T) etc. 

There is, however, an additional feature, impe^rant near 

l £ ^l - A « This is the fee* cf temperature dependent 

occupation of the bare state a d since 

Z = Z • Z 

0 Kcnd o 

where z 0 is the product, cf conduction electron partition 
function and a local level at with correlation U • At 

temperatures * e d , the level will begin to change occupancy, 
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and for T»|s ( n 1 if on-"* the nonmagnetic level will be equally 

2 * . 

occupier* , sc that new, p eff * 3 p 1qwT .If IsJ * A -v t k , then 
this thermal emptying car. lead to serf cue changes in X(T) 
behavicur ever th- r temper-' 'u.-v cctle t k ^ A with Tespect to 
the pure Kcr^c dependence, It p<- rrteular, it is easily sr^n 
that for +ve art. cf order A , x(T) varies non- 

monotonicslly with temperature. This has been found also by 
EWW* They fin' 5 * a peak in x sorting in at ab^uv -v - „ , 

We have net d on e d re* at l~h c&lcul&ti-ns this ^gime since we 
dc net expect, cur connection to be numerically accurate ^h^re. 

Section 4 : THE MIXED VALENCE (^NDEHSON MODEL C ^ C) IN PATH 
INTEGRAL METHOD. 

For this case, the formalism runs parallel to the C = 0 
case - we take the Gaussian integrals in the same way, rake 
terms only uptc 2nd order in g , and arrive at en expression 
for free energy in static fields the minima of which give us 
the HF solutions for C ^ 0 Andersen model. The formal express- 
ions for c/^q. and B a , the adiabatic and transient terms., 
are also the same. The main difference is that -:h~ path now 
fluctuates among three values of £ rather than two. Let us 
denote the thr<^e potentials., corresponding tc “■hr ohree HP 
solutions, by £ £ __ and Then a typical path can be 

shown by figure IS. Note the relation bar-worn pot tale acting 
on opposite spin local electrons.. 
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The long range part cf the transient *--orm contributes 


(expression (32) 


l 0 h = 


1 


t .+T t . +T , r / , \ 

.3 O i o dE (t 1 ) 


2tc 4 


In |t . -t I / dT / dm’ 


u d 


t. 


dT 1 


{ 


5 0 (t) 


C a (x)-S a (t') 


1+Sg( X) 

ln 5 ?|rro } 

Following the same tricks as in (34) ?nd (35) and adding the 
contributions from both spins we find *he coefficient- cf 


(l/m 2 )in |t. - t. j *o be, 

1 J 

(tan -1 ?a - ^n" 1 ^Ktarf 1 ^- ten” 1 5^) + Ct.tn" 1 ^,- 

x (tan" 1 ? ct .^ tan” 1 £ d , ) 



) x 


where a, b, c, d = , £_ or (any of these);, and the 

corresponding s> , b* , c ! , d T = £ ^ £ + or respectively* 

This result is obtained when one cf the local electrons (say 
spin a ) changes pot, r^ti e 1 from ? a to 5^ at the flip i and 
from ? c to at j , and the opposite spin, from C a ,tc5 bl 

at i at « to £ d , at j. 

Thus there are four kinds of exponents, (2 - e ) , ... , 

(2 - ) , three t's corresponding to the thr^ kinds cf 

hops (5 + , £+ * € 0 » £„£ e 0 ) » f!nd thr ^ correspond lag 

J *s . One can write scaling equations involving all these 
parameters, following -knd^rson and Yuval f s method, ¥e assume 
that it is possible to vary o»s in such a way that all the 
e . *s (i = 1, ,,, 4) reach the value of unity simultaneously. 
Even then, we are left with three Jj_ *s, or, in the resonant 
level type of models, three kinds of mixing. These three kinds 
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of mixing could b-e token t- denote either (1) mixing between 
the conduction electrons f n-t either of two resonant levels &r,c! 
between "’he t-wc resonant levels themselves* cr (ii) mixing 
be + 'weo,n the ccr^uctl on electro ?s £ *d either of three ^^scnant 
levels. 

Unfortunately , sclu" i ers tc ?uch “wc end “h^-n ^es^nant 
level mod n is tr~ n r * known. The ref sen we ho no* gc * ntc the 
details cf scaling end the p&romr-ters . 



CHAPTER 7 


BEYOND THE HFA - THE QUASI PARTICLE METHOD 

Section 1: INTRODUCTION 

In Chapter 17 we applied tc rhe non d eg -'rate Andersen 
model with C = 0, the path-^n' ^gral method starting from fche 
H-F solutions. The results gave a strong indication that 
as ■* 0, we have just a strcng-coupleo , Kendo system. However 
traditi inal methods of •'■reading a Kendo mo^l, such as pertur- 
bation in coupling constant J, or scaling methods (again in 
powers of J) b^cam^ more and mc^e unreliable as J approaches 
unity. Therefore, the path-integral method Is me* very useful 
for obtaining quantitative results for physical quantities 
as o. Moreover,, as we briefly showed, the C ^ 3 case 

becomes^. -In the method, qui oe intractable. For similar reasons, 
more realistic or complicated local-electron level schemes 
(e.g., degenerate orbitals, Hund’s rule, CF splitting) became 
difficult to treat. It is also not clear hew a lattice of 

impurities can be dealt with within this treatment. 

, t - ■ - 

1 The method; used in ' this -chapter -to consider these problems 
is an interesting, formulation 'Of many $ ody th^cryv-DeDopinicis ? * 
et al* 25 > 26 showed quite generally that the partition function 
Z 1 of? a large interact! n g system can be written £s, the • 
product .of . the *Z of a noninteracting system and that of; a 
system consisting of noninteracting "quasi particles "-with ■' : , 
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discrete real energy levels. Th® energy levels depend upon the 
interact! ons , ana the op? rational part of th” method is a pre- 
scription for calculating these energy levels, basically 
perturbatively , as a function cf interact! n strength. The 
method has been applied tc a variety cf systems, mostly elec- 
trons scattered from localized impunities' 5 , Applying It to 
nondegenerate Anderson model for small values of [e^l , we find 
that it gives good results for magnetic susceptibility, valence 
fraction etc,, as functions cf temperature. Preliminary work 
or. C ^ 0 case, and. on more than one impurity case show premise. 
This method has the advantage of treating Kendo effect properly, 
as well as the two-impurity problem discussed in Chapter II. We 
therefore hope that, when extended, this method car. overcome 
most cf the difficulties encountered in earlier efforts and 
produce quantitative results for one as w-11 as many impunities. 

In this chapter we firs* briefly outline the formalism and 
then its application to the nondegenrrate d-crbit&l Andersen 
model with U = « • We then discuss th® calculation of quasi- 
particle energies, susceptibility, specific hea*, valence frac- 
tion etc. In the last section the question of interacting 
impurities Is taken up. One fs new tn a position tc discuss the 
competition between Kondc effect, ferro- or ertiferrcm&gnetic 
ordering, and ccheren*' or incoherent mixed valence ordering. 
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S potion 2 : FORMALISM 

Spot! or 2(a) j THE GENERAL QUASIPART I.QLE FORMALISM 

Bloch and De Domini cis start from the perturbation express- 
ion for the partition function - 



» £ 2 
l+l Tr / dr . . . / do 1 e 

n=l o o 


-t3H r 


V( t ) . • • F( ) 


1 + l ‘JJr / dr n . . . / e 
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^ - e Ve 

- "H 0 (r -r 1 ) 
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.P ;2. OH 0 _E 0 T n ^-H 0 (x n -T n _ x ) 
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0 1 ... (i) 


In expression (1) , ( - x^_j) ’ h r - iP* 1 vine lnxr>~val 

between two successive V T s during which the system Is In an 
excited state of the Hamiltonian H 0 . L>~t us denote the energy 
of this state, over that of the initial state, by e i • Than 
the exponential operators in (1) become factors cf the kind 
exp -,))• These time intervals can be easily integrated 

ever and we obtain. 


Z_ 

Z, 


i + I l ("D 

n=l all dia» 
of 0 (el) 


MS jA- 
2u i 


/ dz 


a -pz 2n-l 
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( 2 ) 


where N is the statistical factor (occupation probability) 
associated with the initial soate, and M is The product of 


matrix elements cf V, » 

Regrouping n in (2) 4 n*o products of i -Reducible 

Z ~H 

parts - these separated by a denominator z , or Intermediate 
energy zero, or initial sta*e — it can again be shewn that 
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where nj is the statist-ice 1 weight, cf initial state i and 

is the sum of all irreducible diagrams (each of which contribute 

n ~~ — ~— ) with initial state i . Few Lagrange's th^orem^ 
i e “ e i 

can be applied to (3) and w <-> have 


z/z 0 



... (4) 


whore E^ is determined from the self-ccrsi stent equation 


Si = r i CSj) ... (5) 

The energies E^ are termed quasi particle energies. Now we 
describe a scheme for deriving s Ej for the U = 00 Anderson 
mcdel. 


Section 2(b) * DERIVING THE SELF-CONSISTENT QUA SIP ARTICLE 
EQUATION 

The Anderson Hamiltonian, 


H - k l - l 


can be broken into, 


, U 

' da aa T 2 


F n, n, + V T ( ct c , +c t - c. ) 
A aa d~a v ka dcr da ka' 

a 


V r U 

H o = l e kcT n ka + a e drr n drr + 2 a ""d-a 


v 'ka ka % da da 
ka a 


K 


n , 


+ 


and H* = V I (c? - c, + cT c. ) 

v ka da da ka 7 

xL p O 


... ( 6 ) 
... ( 7 ) 


H 0 above has four initial states according as the d-level is 
singly occupied by either spin (denoted by the subscript a ) , 
Is unoccupied (denoted by subscript 0), or occupied by two 
electrons of both spins. This last state Is excluded If U = 


00 
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Fig 13 , IvreclucjMe dUcjqrqm contri buft'ncjj hnd 


<3 


<1 
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since it then acquires infinitely high energy. The d-elec'rcns 
In the three states have energies e d fend 0 respectively. 

The time sequence of the system (explained after eq. ( 1 )) 

can be expressed in terms cf diagrams. In fig. 13, some cf the 

irreducible diagram contributing to r and r are shown. 

o' 0 

In above diagrams, the dotted l 4 n^s denote occupied d- 
elsctrcn states, spin being indicated on the diagram. The solid 
lines denote electrons or holes according as +he arrow is 
pointing forward or backwards. The joining point between a solid 
and a dotted line contributes e factor cf V, The intermediate 
state energy between two such successive points is determined 
as follows — Draw a vertical line between the two points. Add 
the energies of all the lines pointing right and subtract 
energies of those pointing left. In the diagrams for , if 

this vertical line cuts a dented line, there is no contribution 
fron this line; if it passes through a gap in the ^ct^e-d line f 
a contribution cf - e dQ . • In the diagrams for , this rule 
should be reversed. Th* statistical factor associated with each 
electron or hole line must be written down, and finally summed. 


Following the above rules, we first write down the contri- 
bution from a few diagrams and then expression for r and r 

o 0 

In some approximation. The first diagram in eq. (8) fig. 13 
will contribute., 
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Similarly, the second diagram contributes 
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If we sum all the diagrams in the first row of (8), w^ got 


D l" f ( £ ka) “ 

- i d£ ^ U + nil { 5 I («.«ko)/(^ d< .-0>“] 


'kcr 


= |/ de kc - W Ks^ka'i •■• (I 0 ) 

with I(e,e k<J ) =Z,/ d =t-a' f (e k - a - > /<e+e k , 0 , - e k<J ) . 

Looking at the seen nd row of diagrams we notice that their effe 
is to dress all the internal d-llnes (thcs- occurring between 
the first and the last dot) in the ,«* we y ee for r„ , but 
that new the initial energy for =scfa of these lines is net e 
but e -+• ^ * Sc these taken into account will modify 

expre-a&i on {10} -as> 


r 0 (a) 
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£ dE ka (1 - f(£ ko )) /f E+e a 
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'kcr 


f ( £ k'cr^ da k* a 5 


£+£ k'c t " £ k0 +r a^- £ v 


1' o' 


•7~T ] 
kcr' 


... (11 


There- remains the third row cf diagrams j r , (8) which we shall 
net- take inre account since the corresponding expressions are 
quite complicated and nc suitable way was found to -eke them 
into account. Keiter and Kimball's expression i s rhs same as 
(11), since they use the same approximation. It must however be 
noted, that these could contribute Important logarithmic terms, 
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as they correspond, 3a the large je^| limit, to seme cf Ncziores. 
Do-Dcminici s diagrams^ f cr the Kordc Hamiltonian. Following 
Toulouse^® it can be easily shown that in this limit what we hav' 
accounted for, is Kendo model diagrams of the typo shown in 
fig. 14(a). And what we have net accounted for, is diagrams of 
the type' shown in fig. 14(b). (Here the double line denotes the 
impurity spin; the do+s*, either J or J a ; the conduct! on 
electron or holes, denoted by single solid lines, have spins 
according as the scattering is J z or J ). 


For r Q we make exactly the same approximation and obtain 


the expression, 


r ° (£) -su 


.h 


de n „ f(e kcf )/[e - 6^ + e 


'kcr 


'dcr 


kcr 


A z 11 
- / de 

71 -D 


f ( e k'a' ) 


lc * or 1 


rjs 


J 


'k f a ' w ka 0 k , a l & ka*' 

Comparing. (11) and (12) we find that the integral expressions 


(12) 


in the denominators are nothing but r 0 and ? a respectively. 
Thus (11) and (12) become, 


r cr > 


A 

TC 


fV 1 ~ f(e kg ) 

i ka e+s d a - E ka- r o ri + e da- £ kJ 


• • . ( 13 ) 


r o Ce ) = 


l Ae 


f(< w 


D fecr e ~ e da" l£ ka~ti^ £ ~ £ da +£ kcP 


9 9 9 X 


14) 


Those equations are the same as these obtained by Bringer and 

Luslfeld 30 . 

To determine quasiparticle energies, we must solve the 



self-consistent aquations 
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E_ = 


J d £ 


1 - f < E kP 




... (15) 


E, 


£ y f , f ^ £ ka^ 

- , £ l:a E7^ d ^ to -r^-V dff +e ka T 


71 


... (16) 


section 2(c) • the quasiparticle energies at different 

TEMPERATURE 

2(C)(1) - ANALYTICAL FORMS FOR QUASIPARTICLE ENERGY • 

The self consistent equations (15) and (16) aro difficult tc 
solve end we must make workable spprcximati ons. With this aim 
in mind we first look at the lowest order approximation- to these 
equations, i.a., 


Ba = 


% 


D 


1 ' f < e ka } 


^cr r £ dcr 


Sv 


de 


kc 


. • • (IT ) 


E, 


= - z / 

11 -D 


u 


f(e ka> 


E, 


£ dc + £ ka 


ds 


k a 


• * • ( 18 ) 


At T = 0, fk is a step function and we can perform these 
integrals right away, and obtain,, in "‘■he absence of a magnetic 

field, 


E, 


= r la ( 

TC 


E 0 + E a 


|/B > 


9 ?9 


(19) 


E „ = 


f 4 11 < |e 0 


e d |/JJ} 


( 20 ) 
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Fl<5. \5~ Graphical Solution /o lowesl order yuasi par/icle eyuafii 
ai zero iemp era fu re . 
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The graphical soluhi ~ns to (19) end (20) are shown 3n fig-. 15 


For e -i ~ ve and I s , I > > a * these sdu'ions are approxi- 
a 1 cl 1 1 

mately as follows:- 

From fig* V. 1, w? have only ' , ne sole , glv" 5 - by 


E 


cr 



leg + (A/n) Ir (U d l/D)! 

1 """" ~ 


( 21 ) 


From fig. 15 we have three solutions corresponding to points 
a, b end c 


E 


o 


2_A 

n 


In 


e d + (2 l /%) 


ra ( |e d | /l)| 


* e d 


is 


£ d + D e‘ 


* » • ( 22 ) 


Out of the thr^e roots (22) , we choose the last, one (c) 

B q = £ d “ D exp (7r£ d /2 A) 

This is because at high temperature crly this root continues 
(the ot-h°r two become imaginary)* Besiks, th<~ Integrals (17) 
and (18) should strictly be i n ^c. whose only 

real branch is the left or.e in fig 15(a) or 15 (b). Thus »c f 
will be the only real root* This root rr-pr^s^nts a physical 
bound s-*-ate. The nature cf this bound state beccmes cl~er from 
explicit ground 3tste wave function of Verms an^ Yefat 33 , 
since their ground state also turns mt to be non-m&gnetic with 
energy e d -i) exp(Ti;£ d /2A) • Actually we find above that at 
T c 0 , £3 + S fl < Eq« But we will show in what follows imme- 

diately that the non-magnetic state is indeed lower in energy 


i 
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if we make better approxina fc ions for and. B 0 

To make approximate "rs further than (19) and (20) we act.e 


that T a (E 0 - e 


+ e, ) can be written as 
a.c key 


T o ) + C E n - e ft ^ + e lr 


do' 


Til 


) (Ecr ) + 


... (23) 


If wo substitute this (16) then for relevant- values of 
£ ka (- ve st T = 0) and for E c an- 1 S a fourth ab-ve* the 
first term in (23) dominates, So we approximate 




= l 

a 


A 

% 


! 

■1 


D 


f(e 


ka ; 


Cl£ 


o 


£ da* 1 ’ 


ka 
£ 


ka 


1L 


(24) 


But a similar replacement cannot be dene in (15) sine* 1 
at E 0 t T q (E) varies quite rapidly* So we leave this equation 
at its lowest approxim&ti Xj as ir. (17) i*e», 



A 

71 


/ 

“D 


1 


.±£fkal 

E a + 6 da 


'ka 


■* * » 


(2s: 


Bring© r snd Lust f eld also started wi*h the same express ions (13) 
and (14) for E 0 and E ^ • But their r, esul~s are poor f or good 

moment (low T^) because of ^his self-energy addition 1. n the 
sensitive denominator. 


These approximations to quasi particle equations we^e 
derived starting fr^m the solutions their lowest order 
approximation at T = 0, At high (high to be clarified later) 
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temperatures, we can make another kind of simplification. For 
this we truncate the quasipert Icle equation (11) at 


E fl = | (l-f(£ v „))/[B 0 + 


A 

TC 


ko ^" MC k^ //Li3 a 1 °dcr 
I / de ] r j ^ t f ( £ ^ ? ^ i ) / ( e 


'kcr 


« « » (26) 


a 1 


a k y a 


*' C 


key 


) 3 


This corresponds tc diagrams of the class shewn in fig. 16(a). 

Toulouse^ has shown that for large j s I this corresponds to 

a 

the diagrams cf fig. 16(b) in the Kendo model. In 16(b) double 

23 

line represents the K n ndo spin. Noziere’s and D --Domini els 
have shewn for the Kendo Model that these arc the diagrams 
contributing most divergent logrithmic terms. So the approximate 
form (26) can be said to be accurate upto , *.h?s order. 


We similarly truncate (12) at 


= I i / de lm f (0/!> 


kcr 


'kcr 


a + ai- 
der 1 


G 


~ n f a£ k 1 a ’ ( 1 " f (e k i cy «))/(li a “£ k » a .+£ kc y)] 


. .. (27) 


We now calculate E c , B a etc. from the equstl ons (24, 25, 26, 
27) vhich Involve ‘rregrels cf the form, 


1(a) 


k f(x) dx 

T \ x + a 


This integral has been evaluated in appendix VI and we quote 
the result here 


IO) = In --fM- + K [In (m/2ipa) + 2^ (ipa/m) - i|» ( ipa/2n;) ] 


1 
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The low end high temperature expansions of "this are as follows:- 
1(a) = In -L~i- - (™) 2 » it! << a 

1(a) 1 In (q |) + l^( 3 )(|~) 2 y 7tl» a 

where q = ~ e^” r ^ » 


Sec. 2(c) (11) i LOW TEMPERATURE LIMIT 


We now leek at quasiparticle energies 1 o the low tempera- 
ture limit, first B a and then E 0 - Equations (24) and (25) 
arq. 


% = 


_ 4 

u 


! 

-D 


D f ( a ) ds 

iTF d rE 7 


A 

TC 


1 < e d +:B o ) 


• « • (24 ) 





D 

/ 




f ( e ) ds 


e+B, 


■ £ d“ E cy 


2A 

% 




• % • (2 5 3 


In the absence of a magnetic field* 

At zero temperature wo know those integrals rc be of logarithmic 


f erm , i*e. , 


E 


a 




v 


E 


o 




The solutions to those are, 

•a l £ d + | ln UalAW 

E -V f in — i± a 

a — it D 


• • • (26) 


Eo i.e d + E a - D exp (xE Q /2A) % e & -^ a -fDj^\ exp(Tce d /2A) (27) 



103 


i in (27) is the 

The reason for choosing only the root m--n 

. , -u r^SeiPt f* 1 cl« SH'-TgleS 

seme as ?xpla*n*d in sec. 2(c)(1). Th- Q u 

«• JS 


are, E 0 


and E^ + B ( 


(= f ) f0T .noTmcgne-*fc and 


magnetic state? respectively. So from (27) v 


t = ~ f c 


i~T exp (Tte a /2 A > 


(28) 


. . , , r^s indeed lie below 

Thus the nonmagnetic quasipar^icle leva -l 

A ^ , T , -nail later identify 
th* magnetic cn?, by an emeunt that v- 

with T k , the Kondo temperature, 

, „ , -c- and in particular 

Now, in order to see how E 0 cino 

-.rnopratursS'. we differ— 
the gap , change at finite but lew *•>* 

entiate equations (24) and (25) - 


A d_ 
ir ■ d.1 


f f ( e) de 

f J e + e, + E 


II ( i‘ a i‘^ 2 , 2/[i- 1 A(ea +B a > J 


S 2 = I 5 fi^ dE f(0/(e+E 0 -e a -E o ) 2. - r \ 
nnglecting a/ e ,= i' k /4 In comparison wit.b ( >» 

tt -- VP’ have used the 

is the limit cf zero temperature. Here - 


% 2 T 


(28) 


*;• (29) 


3 E 


standard relations-, 


L f i§I 

31 


3x(e) 


T 3e 


ana /“0 (e) |f de = -0(0) - {£ T 2 f-f UJ 1 ’* 0 

— <o O £ 


(30) 


( 31 ) 
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Same results are obtained free: ah' low temper attire expansion 
of 1(a) given in appendix VI. 


From (28) end (29) we see that as temperature starts 
increasing from zero.-, ih? m£gni > ude of the ^-r Tgy gap 
S D - B a ( = H ) Increases. But as has bee-; show’- graphically 
In fig. 17 after some oemp^ratu^e It deceases, becomes z<=ro, 
and ultimately becomes positive, i.e., nonm&gni^u^e quasi - 
particle state lies above * he magnetic quasi parti cl? state. Tc 
see this, tsk® up the equation 

E q = ;j:»* I (S 0 - E a ) 

The graphical solutions to this are shewn in fig. 17 for 
different temperatures.- 

The curves are plots of - 2 --- I (3„ - f j , T-, < T 0 

% 0 7 ± & 


From this we c&n also obtain oh^ value of T at which the gap 

* 

closes^ i.e., the n* nmagne-J c and the magnetic quest pa rticles 
have the same energy. Th'S will happen when the straight line 


in fig. 17 cuts the curve a--, E^ , i.e., 



2a 


7C 


KO) 


But w e knew th£t 1(0) = In (q.^) (sc-" 5 high +'--rap' fc ure express- 

ion cf 1(a)). The required temperature will thus be given by. 


= 


or 


T = 


2 A 

it 


111 


til 


1 

1 


1 ext 


M v 


:a )/24] i|vrf Sd > 




: d /2A 


• 


(32) 
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since S does not change much over this temperature range 

C -v T k ). 

Thus we have shown ^hf-t st zero temp- mature, the ground 

ste+; e is that of the ncnmegn=ti c quasi particle which lies below 

the magnetic quasi particle levels by an amount Tg • This state 

is a mixture of nat = 0 fend equal amounts c f -he two n<f = 1 

states with mere cf the latter, as is made clear fr"m the valence 

fraction obtained (sec. 4)* The difference between rh» magnetic 

a nd nonmagnetic qu&siparticle levels vanishes ct a temp ^ , 

q 

a t. temperature above which, the* magnetic state rfo r ground 
state* At temperature much higher than ‘“his, as we will show 
below, this difference is OCje^} •)• 


Section 2(c) (iii) ~ QUASIPARTICLE ENERGIES AT HIGH TEMPERATURES 
We now turn to Ob' 3 high temperature equations (26) and 

(27) . 


If for ns a 
term cf high 


-D e k + S a+Ea - fl 
- e k ) ;n -hr 
t empers tu^e e X pars 1 on , 


... (26) 

<W 

wr simply put the first- 
then we ge + , 


* ~ I (E + £ , 
a — % K a a 


f- m ( 4 $)) 


(33) 


egajn applying the. high temperature expansion, 


” m (fr } + } 5(3) (-^ 


. 2A qT 

+£-, - ~ — In ■?“— 


'd n 


1 


tt;T 


% 


In 


(f) 


+ } 5(3) t(e a - A 


In $~l /itT)) 2 


(34) 
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Similarly, fcr S 0 , from equati :n (27), 


1 ” m J c(3) {( 


% 




.. (35) 


Th^ assumptions going into (33), (34) and (35) ar~, 
I » E a , T» P in (q ^-) 

These are satisfied for T » a 
F rom above relations we find that 



This shows that the difference between the magnetic and non- 
magnetic quasiparticle levels starts diminishing as ofte tempera^ 
ture is lowered* This is consistent with the low temperature 
results. 


Section 3. THE MAGNETIC SUSCEPTIBILITY 

We have seen that ?he partition function of the system 
resembles that of a system with three discr^e levels with 
magnetic moments ccrrespondi ng tc + 1 or ^-spir. spates. There- 
fore we expect a Van-Vleck like susceptibility wh^re the levels 
shift with temperature — one in which low temperature suscep- 
tibility sh-iws zero moment (that of the ground state) and high 
temperature susceptibility shews local mcmen < *. Moreover, the 
magnetic moments of the quasipart. teles , or spinlike states do 
not remain constant and vary with temperature. Therefore we 
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obtain th~ static susceptibility mme c ly by censi-- ng '■‘Yip 
system in e magnetic fml- 3 -. Th« pert in or fun erf •*> is-, 


2 = | esp [“P(e da +b a )J -i- exp (-po Q ) 


••• (36) 


wh-vp e 


da 


*d 


P' d c>H 


quasi pert icle energies "i 
The static susceptibi liny 


, and £ - "i'Z; and E are the 
" cl a a o 

the p’-'SPnco of a ms gn^ic field H. 
i-s giv-^n by -'h -7 c-xp^-ssi'T , 


x 


1 9 2 In Z 

P 3 1-I 2 


Per rh'*- partition function cf (36)-, w? ge*-, 

-g(e d -B 0 ) 


X = 2{ p(op d -B;) 2 -B")^ 


,QX‘ 


. lT1 t; e 
" Z 


... (37) 


whore primes and double primes dene be first and second 
derivatives, respectively, w»rv. L v H* In (37) w-= b&v-^ put 
E 0 f = 0 which will be shewn b e lew * 


Se<? tint* 3(a) i 2350- LEM PEE A'f IG3 SUSCEPTIBILITY 

To find the field -’ariva-ives r f Eq an -1 s ff we again 
take up their lew end high ?mp?mur<. =xp as si “ns separately. 
From the. lew temperature quasi parti cle squall ms (24) and 
(25) we have 


E = 


A 

% 


-l 


{1 - f(e ka ))ae ko 


A p d-f(e k ))fc k 


■^cr" 8 darker 


% 


I 


-L ^cr'^d* 


■^ d H- 




I (» 


'd 




4*4 C3B) 
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E 


„=•£!/ 

'TT * 


^ ^^ e lca^ de k 


71 o -X. B o- e do +E kff- Ji a 


A. / 

71 J 


V f( £-,..) dj 


k 


T , E-e.+a^H+e,^ 

1> 0 Cl Q. lv CT 


= ~ I I(E "E -e-,+aa-,Pl) 
jl O 0 G. a 

<3 


... (39) 


whe^e we have put e, = e. 


crp c H transformed the 


'kcr " “k 

i’^-'gretion variable i o s-^ which intro^uc^s H-dopendonce at 
the banc 1 edges which can b^ negl»c*e--i • p c and \i^ are the 
magnetic moments of birr conduction end c* -.electrons respec- 
tively. 


Differentiating above equations for and E 0 , we find, 


K =n ( K - °l J d) l'(V e d-^d® 


(40) 


E- 


A 

7C 


J (B ' 0 - s; + 


... (41) 


where I* (a) = da 1 (a) 

From (40) end (41) we find that as H + 0 


E 


, *1 (%)/* 

= _ cyu - - 

U 1 - Al’ (E ) /% 


and E 0 = 0 

where Ecr = + E ff 

Differentiating (40) and (41) tree more, we go 1 


... (42) 


E 


a 


% 


(E;- atld ) 2 i" ay + 4 - b" i' ay 


(43) 


e: 


I (E 0 -E c7 ) + (E”-E")I i } ... ( 44) 


tc - — 


o cr 
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Substituting from (42) int? th^se t nd rearranging, we find, 
2 A 


■n>" 


= '~d % 


I M (js 0 ) 


{ 1 


71 <7 


9 


.. (45) 


B'A = 


P.| 4 /it 


(i -§^’(-o-2„)hi-^i : (S„) J 


'0 <7 


A I ! (E ~B ) I"(E ) 

5f [2I"CV*W> ° ° ° 


T t.. (46) 


Substituting (42), (44) an* (45) in (37), th« susceptibility 
can be found. 


We first find x et zero temp-rc-ture. For large l e ^l 
zero temperature we found that 



, v £ fi + 
— Cl 




D 


at 


and E 0 - E^ = ™ 'fD Je^J exp (Tie d _/2A) 

Also, at T = 0, 1(e) 2l ln (|a|/D) 
Therefore at this temperature, 

If (a) = i and I" (a)' = - jk 

c* ci " 


Since at T = 0 only *-h« nonmagnetic stfte is occupied, from 
(37) and (46) we have 


X = - Ee" = 


» a d 

(B 0 - ® c > 


4 .. (47) 


where we have neglected terms o(A/7t|e d | ) and 0(i' v / A ) in 


k / 


comparison with unity. Thus we have, 

2 

^d 

X (T = 0) = 


(47 f ) 


VD |e^j exp(7ie d /2 A) 

We can now identify the denominator in (47 r ) with the Kendo 
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temperature since the result (47’) is standard end is, in 

form and upto numerical factors, obtained by all methods appli- 
cable tc the Kondo problem at low rpmperatures — scaling applied 
tc the asymptotic (long tine) approximation*^ , numerical calcula- 
tions applying renormalize r.i on group rc^thc^s^, *-tc. Thus we have, 

T k = V'Dje^'l exp (me d /'2A ) ... (48) 

= “ exp (-1/JP) 

T o 

(This larger is the standard Kc-rfo expression). 

Thus, for the good moment- ( |e d | > > A ) case we identify 
the effective exchange coupling, J, as, 

r 2 A 

1 I l e aT 

This is the same as obtained by KWW. 

Haldare has pointed out that itfcx* t 0 should be IFAD” f° r 

a U '« co model rather than V D [e^j which is obtained by us. Thus 

our expression is deficient by a factor ] [~r^~T v \ r I j Tp * ^tmytage 34 

rdl 

has shown in the scaling approach for th= Kon^ o model that- the 
contribution comes about from III order scatterings of the type 
shown in fig. 18(a), which will contribute to Z diagrams of 
the type shown In fig. 18(b). In the Anderson model, in our for- 
malism, this corresponds to the diagram in fig. 18(c). As men- 
tioned in sec. 2, wp have not takers .such processes i.ot-o account, 
and henc-^ the missing factor of YJ"J | p . 
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The abev- oxpfessi have h^'t 


, r 4 V d rely ff-T |e I » A • 


i. n« nc W® 1 !oV0 

But- th^ identified or- with Ik ^ 11s u “' ,h ~ ' 

j. . ^-v~ran£ (except these 

taken Into account r br meat imp j -■'•■••■ •- 

, v4 r f - 1 ? rams should give 

already mentioned). So the Sf.m j£. 

„ i i > •. , -i^u u^h in this c?so the 

correct results for lap 6 ^ » 1 51 J J ' 

i= /•-■’i energies 

expressions are no lorg-r so c~npcC. ( c - 

a tT _ ^ r ' r * rx S'"‘ ft tO 

bPing rf hh’ tw order, ' 0(M)t w hav ' ■ 
num prl c& 1 c& leu Is ti -hb • 

S^P 3(b) : THE MAGNETIC 3U3CEPT IBILITf J HI® TEMPERS 

TO cbT.alo an gross' ~: for x .* Wgh WW we 

4 , a, -hi -*Yi fh-ro^rs^u^ quesipt rti clc 0QUc> 1 
l^k^wlse with the high re jp-.rc ,.u ^ 

Before that, v, »r «*•«*' ou-iparti ele level Is 

mu ch above ( * W 4 I > «»" “«»««> 1=Y ° 1S ’ ** 

, .. i-vio-ppf nrp contribution 

J t -,i > > 4 . S° its occupation, *nP *-n- , 

to Z or X can bo neglec+c-d and w- may wrt ~ r 


Z 1 1 


C49) 


U 

j, , eH^T^y which 

so we need to know only *he magnetic quesip-- - c 

is given at high temperatures by 


% = £ / 
% \ 


(1 - f v ) da,. 


^ Id, 


: d - 0 p a H- 


I — v™— — -or - % a a • u. n u 

li E 0 +s d -op- d H -- r~'l n ( lj)) "' £ ]r ••« (50) 


whew we have -ltalnatert P 0 -lb the same way as for aero 


t em pe rature c a s ; 
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Prom (49) we obtain the expression for 

-pi 


x as 


X = 


| I (BJ - EKc^-E;) 2 ) e 


CT 


Tjtf t j- n 




■an ' 
a 


... (51) 


Differentiating (50) twice, in the limit- H = 0, we obtain, 


E * = 


- I' (1* 
%■ cr 


2 A q.T> 


7t 


a In 


Al'(l -24-in P) 

7C ' C I JJ ' 


... (52) 


2nd 


O A I ,! (E rt -|~ In 4^ 
„ 2 r A cr - it jj 

E " = C 


a ' a 7t U~~r (I a -.pln } 2 


TC 


An 2 


I, tVr Ul ^"(iP-Aln B 


+ ^ (1 - P' (E In P) 

7 C rf it. t 'l J s 


m • • (53) 


At sufficiently high temperature, but T < < I e d | » (50) can be 
written as , 

. I £ n + E - -- In 1 

A . 1 a a 7c i , I 

D 


E o =“ In 


% 


So that 

E a * 1 C 1 + 


£ . + i m ~ 2 a 5 t 

CL -,o " n D 71 ln D 


and 


A f,i 


E " 2L “IJ-d ]7 I 

■ r , A l°d 1 2A _ qlN 

( £ ^ 1 tr In . n • - 07" In rp) 


l\2 


'd ‘ tc D 
neglecting terms o(A/|e d J). 


% 
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Punting ^hsss in (51) g~t the result, 


xi Puf [1 + 


2 A/ue - 


a 


1 + — — In 
TC£ d 


TUT" 


D 


2 A 

ire , D 
Cl 


This is the standard high temperature result for a Kendo system 
with J p replace by the value 2h/%e a , The expression (54) blows 


up st 


where Tg 



= VdT 


*e d /2A 


Cl 


We h&ve 1 / q = .882 


this factor should be .303* 


identified in the previous section, 
whereas Nzier^s has pointed cut that 


3(c) . - NUMERICAL RESULTS 

life have discuss'd above analytically the susceptibility for 
high and low temperature, especially for the g-oc* moment case 
(a d -ve , je d |>>A). In addition we find X (T=0) for all e d . 

In the general case, I*e.-, fer all temperatures and e d i x(T) 
has been calculated numerically by solving self consistently for 
the quasi particle energies E 0 and E a (separately at each 
temperature and for each value cf a a ) , and obtaining their 
field derivatives. We exhibit the results cf seme numerical 
computations in fig* 20. 

On comparing these with KWW we find that cur susceptibility 
results ere qualitatively correct, but higher by a factor of 
about 4* Also, nenme net on! city is more pronounce** in the e d = 0 
case. Our model is also somewhat different from KWW as we have 
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taken U — ~ whereas ■’•hey have r sken U = D/2. Si. ice tii'- 
ch&rectcri sti cs cf these cunr^s l re very si ml If r 4 c KWW, wa 
expect similar deviatl :n? from universal! ty. 


Section 4 : THE VALENCE FRICTION 

The d-electrcn cecupa t-1 on, tt , as w? call It, ^hr valence 
fraction is a quantity which h^lps us In understanding the nature 
of the state. It is defin'd as, 

■ ' 7=1 <n ia > 

Cf 

Since Z = Tr { exp (-p[e d I * dc + 1 ^ r A/(\ f i c t/r!/ h «c.) 


a 


ka 


:a ka v "led ka^da 


we have 


I <n 




1 8 3F 


+ ^d-Ad+OJ 


d *0' =-?*-(**) =ae a 
where F Is the free energy. 


At zero temperature, F= E 0 , the energy cf the nemnt gnetic ground 
st at©. 

3E„ 


V(T = 0) = 


3e - 
a 


From low temperature equ<-*i "n? ( 24 ) and (2?) we h&ve, 


E o = ~ ^ E o - £ d - E a ) 

®a = iT" X CS a " f £ d^ 

Differentiating these w.r.t. , we find 


!£d = aI ’ (j V £ a)/* 

5e d 1-61' (E a +C a ) A 
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TCl\ 


3E n 3E rt 2AI 1 (S --£ i--E ) /% 

2. = „( 1+ _^L) ^ (l -) ( i 

3£ d 3e d 1-2 AI* (E 0 -£ d -B a ) A 


A v 


m 


-) 


* * * 058) 


'Cl 


putting I f (&) = ^ , E 0 - = - T k , anr 5 taking lowest order 

in 


k A 


, T^T‘ 


The expression (58) 4.? very close to uT*ty. W-> therefore 
conclude that the ground state almost f-nt 1 ly consists of d~ 
electrons. This is confirmed by th° ground-state calculation of 
Ysfet £.nd Varma', Also notice that naive HP A would give 


V * 1 



a2 


) T 


Thus the Ionic model is close to truth and 

very small, * 0 (— •) . 

A 


the relative error is 


Section 5 : SUGGESTIONS ON INCLUSION OF SCREENING TERMS 


AND OF INTERACTIONS BETWEEN IONS 

The inclusion of a screening torn, g(n df +n d+ ) £ c 

kk ! a 

gives rise to additional diagrams of ~he typ* shown in fig. 2 la, 
And all combinations of such conduction electron parts Inserted 
in all the d-lines in all diagrams for 2> • These are self 


energy insertions and have the effect of simply adding 9 

in second order, 

2 r f k ^ " f k'^ de k d£ k’ 

Pg J — ■ 


£ k £ k* 


This has been shewn to be C. In higher order the contributions are 
more complicated and provide logarithmic terms. 
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( 4 > 


diagram 1"wo non- mlerac'Hn^j imp. 

Ft ^ 2.1 
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The other class cf diagrams involving g are cf the type 
that are difficult to handle for the same reasons as mentioned 
in section 2. These diagrams are shown ir fig. 2lfc. 

Thus, right now, we are not in a position rc fully discuss 
the effect of a screening term. However, taking inly the second 
order term shown above, we got the addition of C tc which 

is physically plausible. 

The Interaction between two impurities, with an added 
hopping term, X c &2a + ■ h -° c ” ? precede vie two mechanisms - 

the direct hopping and that via conduction electrons. These will 
give rise to diagrams cf the type shown in fig. 21c, where the 
wavy lines denote direct hopping, X • 

iis for HFA, the quasiparticl^s here will be cf the type (00 > , 
|0c > , | ac> end (c-cr > , but their nature, obviously, is very 

much different from HP ores, since 0 here denotes the ground 
state and is net like the HF nonmagnetic state (which arises only 
for C ^ 0 case). However, once we obtain the Interaction energies 
in the quasiparticle formalism for these product states, the 
pseadospin analogy can bo carried ever. 

We do net proceed with the formalism for interacting impuri- 
ties since, in this formalism, we find It difficult to show the 
basic relation for two nentoff eraettng impurities (diagrams shown 

fig. 21d) 
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% = 2B 1 

where Eg artf ^ are ro^al quasi parti cl-' ^.r-rgies for two an<? 
one impurities respectively.. 



CHAPTER VI 


CCiTCLUSl OHS 


The single valence fluctuation i “‘.purity has been treated 
in the Ear tree -lock approximation, the path integral, and the 
quasiparticle fonnalism. They give successively sophisticated 
results. Only in the first of these (HP A) we have taken 
interaction of the impuriry with conduction electron or 
boson fields into account. ¥e obtain no new qualitative 
results in the Eh' A, out cur approximate analytical results are 
shorn to be quite accurate $ so that theories, which have IIP A 
as basis (e.g. those in chapters III and IV) can get the HP 
energies quickly without going into detailed solutions of HP 
equations for a given set of parameters. 


The path-integral method establishes firmly connections 
of the Anderson model with the Hondo model. But . as is well 
known (e.g. from HW calculations showing deviations from 
universality) , as the Anderson model parameters approach the 
nonmagnetic regime, at higher temperatures, the correspondence 
with Hondo model is not complete. The partition function 
obtained from the path integral method has to be scaled to 
a known model in order to obtain quantitative results for 
physical properties. Scaling becomes tedious at higher 
temperatures (not discusses, in the thesis). 
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In the c unsiparticle metliou of chapter V we believe to 
have found a ns tied which will deal succes scully with all 
temperature and parameter ranges. iho results in this chapter 
have matched correct^ with known results tor relevant para- 
meter ranges. It however remains to treat the model with 
realistic moaifications. 

Iho concentrated IiY system lias boon treated only in the 
HFA. In the no an ficla approximation ic shows an interesting 
new kina of phase, the coherent valence fluctuation phase, 
because of the limitations of 1 -IlTA we nave not cesrriea the 
treatment in chapter III beyond the mean field approximation. 

It however shows the pseudospin models to be realistic ana easy 
to handle . 

'Ihus we find the quasiparticle method of chapter Y to be 
the most useful for solving the single LV ion problem. Once 
this has been done and the two-ion problem has been sorted out, 
the results can form the basis of a pseudospin model for the 
concentrated system, as in chapter III. 



APPENDIX I 


FREE ENERGY FOR 1H3 AIT DIRS OH HOFEI- (C 0) IN HFA AND 
CONDITION TOR KINIhlM ENERGY SOLUTIONS 


H = 7 e, cf c, + 7 £jC + c + V 7 (k c . + 1 . c. ) + 

L ° Irrr Vrr 1<*rr L d nrr G.nr f v ko CLO 

ko 


r - 'ko^ka^ko 1 f °d do* ao 

ICO o 


+ u n d rn “ c(n d+ + n d+ ) / 2 (AI.l) 

'The HFA consists of replacing a pain of operators AS as 

AB ■* A<S> + B<A> - <A><B> (AI.2) 

Doing this to the last two terms in (AI.l) we get 


o 

eff , 0 


H S" + I (< J1 dt >+<n a+ >) 


0<n d+ ><n a+ > 


(AI.3) 


where H^ f = I 4o c da + |/ka °Ka°k a + 7 I (c ka°ka +11 -°- } 

0 Ea to (A.I.4) 


(U-0) <n d _ a >- 0<n dn > . 


where e” 1 = e i(J + 

*f 

Denote s dcJ . = x^ where a = ± 
and < n d0 > = n^ and write (AI. 4) as 
*e = + - C n a 


a a 


( AI . 4' ) 


_p_p 

Kt t is the Hamiltonian for two d electrons of opposite spin 
H-x 1 ' 

with effective d state energy x + and hybridizing with the 
conduction electrons via the last term in (AI.3). The Green's 
function for them will be 


<W“> = (“> + l (“) 7 <\- n (“) 7 W 


do 


1: 


dO 


ko 


'do' 


a) 




l-Y 2 I a k<J (e) 

(ai;5) 


where G d(7 (*«>) 


1 'iff 3116 G ko (w) = 


(D - 8, 


'ko 


(0 — 
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In finding the sum I G, ( w ) T we make the- usual approximation 
"" no 


k 


of summing the • 5 - function part dene • nd obtain £ Gkcr ( u ) 


iA , with A = 7t p V . Than, 
Gd (“ ) - — 


10 — 


1 

feff 

'dcr 


i A 


The density of states for such a propagator is 


p ao (0)) te 


___ A/jl , n , 

- — ‘ £' ' ”* '2 

(“ “ e da > + 4 


and the energy is given at zero temperature by 


/ p „ (& ) • 6J 

' ac v 


-D 


ci to 


lc 


UUS) 


(11,7) 


To obtain free energy of the system we must adc energies for 
electrons of both spin and to it the remaining terras in (AI.2), 
So finally, 


F W = 1 


4 I / 


a -£) ( a) ^ x ) -I- ^ 


a) d a) TT , C t , \ 2 

“* „ jj 2i jx + -at ( n , r n ) 

2 j. A 2 " r 2 T- “ 


a 


"H „ 


UI.8) 


_A f° __ to dt o 

7t ti / \ 2 1 1 2 

-D ( x_^ ) ~ r & 


+ 


-t / df* z-ii ! ( s) 
-S-x, 


-U \ 7 " V* 

r U 


(*J 


= - / aZ-z-a' (s)+ x_. n(x_,_) 




where n * ( z ) =* 


D+x, 


- A/tc 

~~2\ 2’ 
Z + A 


The free energy expression <11. 8) ccn then be written in the 
form 
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%F =“ I * / dZ-z-n'(z) - x n(x )) - U n(x_, ) n(x_) + 
a 

+ ~ {n(x + ) + n(x__)}^ where a - + (hl*10) 

In (;j.8) we had used the abbreviation 

2 + = e d + (U-C) n(x_) - C jb {:-:*) 

Putting xj. back in terras of n(;-: — ) in (AJ.10) 

A 

*HF = - l * / dZ-z-n'(z) - e d n(x a )} + U n(x.,J n( J -- 

CC oo 

- ~ (n(x + ) + n(x_)} ^ t a = + (hl*ll) 

This is the free energy expression obtained by Haldane. 

Now differentiate (~1 . 11} w.r.t* x+ 

9 F -j-j-p 

= „ x+ n' (x_j_) + U n ' ( x + ) n(x T ) - C n : ( x + ) (n(x + )+n( xj ) + 

+ e d n'(x + ) = n ! (x + ) [e d -x + +(U-G)n(x„) -C n(x + )] 

UI . 12 ) 

Sxtrema of the free energy function occur when (hi ,12) vanish i.e#, 
Xj. = e d + (U-C} n(x ) - C n(3^> 

These sre the same as -®elf consistency conditions obtained in 
aquations (*) and C») in Chapter II* 

To find out the condition for these extrema to be minima, 
we must find the matrix of second, derivatives of the 

determinant of which must be + VG for solutions which are absolute 
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minima, i* e* , stall 3 solutions* From (,11.12) 


2 

rf 


-jr- = n” (x + ) [e cL -x + +(!'a-G)n(x^)-On(x + ) ]+n 5 (x + ) [-1-cn 1 (x + )] 


3 X 


'+ 


= -n' (x ) (1 + on' (x, ) ) 

t *r 


(AI.13) 


at an extremum since there the quantity within the first bracket 
vanishes* also from (hi *12, 


3 2 I- 


HP 


ax X = (U“C) n*(xj n ’ ( x + ) 

4 - — 


The determinant , t hen , is, 


UI*14) 


n' ( x + ) n. ' (xj (1+Cn' (x + ))(l+Cn'(xJ).-(U-C) 2 {n* (xjn* ( x + ) Hal * 15 ) 

Since »Hx+)* n (x_.) are both - V3 , the sign of the determinant 
is determined by the sign o'f 

U + C n* <x + )) (1 + C n*(0) - (U-C) 2 n»(x-J n«(i + ) 

= 1 + C(n f <x + ) + n * { ::_) } - '{(U-C) 2 }- C 2 n*(x + ) nHxJ 

(11,15) 

The other condition for minimum is that the trace of the 
matrix of second, derivatives must be positive, which gives us, 
from (11.13,, 

-n»(x + ) (1 + G n»(x + )) - n'(0 (1+0 n« (xj } (11,16) 

Thus the conditions for minimum are that the expressions in 
(Al,15) and (11,16) both are positive* 



APP3HDI2 II 


10 SHOW THAT THE RAMOS OS UWoTAELS SOLUTIONS TO HF EQUATIONS IS 
BOUNDED BY EXTREME III £_ . ' • ' 


The HF equations become just stable (or unstable) as the 
determinant of second derivatives becomes just equal to zero* 


From equation (2’) of chapter II, 


£ d ~ 3+ “ CH - C) n(:^ ) + C n(x.) 


(All, 1) 


The derivative of e d w.r.t. x+ is ? 
de 

~ = 1 - (IJ - C) n»00 , 


ta 4 = 1 - 01 - c) m<o + C mco 


From equation (5») of chapter li s 

2^ = . E d + (U ~ Cv n(x + ) - C n(;0 

differentiating- which, vr-r get 


cU, . 

dx+ 


ae d 

4st— 


+ (U - C) n*(3+) 


1 + G 


(All. 2) 


(All. 3) 


Substituting this in (All. 2) u.i get * 

de d 

Tiu — + " u 0*4x4.) - = :( l'+..G nU 2 +)} {,1 + C n l {::J > 

- (U-C) 2 m(0 m(0 • OJI.4) 

4 V * 

The r*h*s* in (All .4) us the same as the determinant of second 


« »< *> *. 


de 


derivatives of- F^-. Thus 1 when the f *h*s; vanishes, Tnry ■ must 





vanish since it does not guarantee that 
brackets on the l.h, s# does. 
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the quantity within 


11 so . 



de , dx__ 
d x~" clY 


From the form of Y in ecj. (-*1.4'} we see that 
monotonically with x + or x_ . Thus 


it varies 


de 


d 


dx 



= 0 


0 ==> 



APPENDIX III' 

ENERGY OF WO INTERACTING IONS IN HFA 

The general expression for the energy of ‘two ions is 
simply the sura of that for each one of them, in HFA - 


F = l — Jw-dw G? («) - I 
i , o i 


— 4 s 4- 

U ir. n: 
1 x 





<AIII*1) 


where i denotes the ion site so t hut i = 1, 2 , end cl (oj) 
is to be determined for each interacting ion in the HFA. G? (<o ) 
can be determined in a straightforward manner in terms of 

» the propagation expressions for noni.nteracting ions 
in HFA. Doing this and using the HF self consistent equations 
mentioned in Chapter II, we get. 


F = 


A 

TC 


I [ /.— 

a 


ede 


D { e-e,,-(U-C)n a +cii a +y> 2 +A 2 


+ / 


o 


ede 


D {e-e d -(U-C)n +cn -V> +A‘ 


2(U~C)(n + n + + d + <£) 


+ § (n? + d? + n+ + df) 


(.till *2) 
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r* 4- 

l Zlr C l-~ C 
k 


+ 1 


k u k0 u k0 T k e d : n d + TJn ct+ n d+“ 2 Ka &^ n d* 


(n M + n,.)' 


The sum of last two terms cor be written &s (rider *s are still the 
eigenvalues) , 

~T“ ^ n d+ +n d4-^ " f ^ n d+” n d+^ 2 
Now th<= corresponding exponential factor's can be replaced by 
Gaussian integrals in the same way as in the text sc that we 
have d electrons of spin cr travelling 4 n potent :el iY + oX 
and Gausian factors with exponents ._f the form 


-p cf 


~ 1 
U-2C ' 


Sc the total energy of the system, when X and Y ere time 
independent, is the sum of resonant energies of the two d 
■electrons plus contributions from the Gaussian factor s 


A 

rc 




ede 


v-2 

+ tr + tji"2c 


Lilli 


(AIV.3) 


c ~D (e-e d -iY- 0 X) 2 +A 2 
differentiation the above expression for F w.r».t. X and Y 
respectively, and equating to z^ro to fir* its extrema-, we get 

AU 


2X 


71 


[( 


n -l e d +iY " X v it i e.+iY-J-X 

f -tan" 1 )-($ -tan- 1 - a - )1 

^ A ^ A 


(AIV.4) 


p-t y _ A(U“2C) r / tx . -1 e d +YY ” X x , ir -1 

21 Y - — [ (p -tan + (w -tan - 1 - 


e d tiY~X 


■5~) ] CAIV.5) 


Adding and subtracting (AIV.4) and (AIV-.5), and putting 
x+ = iY + X-, 

we get equations identical to the HF self consistent ones - 

x + = (U - C) nCx^. ) - C n(x + ) 



APPEHjjIX V 


THE INTEGRALS OC CURING I IT THE SHORT RAH 03 THRL IN TEE 
PATH INTEGRAL METHDh ALT! EH TEiu : I2TATI 027 CP r 


The short range term, obtained in the tent is, 


1 


Hie 

short 

t . 

t . -rT 

X 0 

X 

/ dx 

J &T 

t . 

t . 

X 

X 


dC ( t ' ) d 


1 ? 


- [-jr— £4* initCLd-] 

€ 2 (t)-^(t') 1+5 2 (t j ) 

(v.l) 


Integrating over t, if £ + ano are values ox £ at t_. -i-x and 


t. respectively; 


i o 


S = / dr* - d -4fP- { in i V T o" T ’ I £~ 


1+5 


+ 


sf-? 2 ^') 1 1+5 2 (t ! ) 


In 1 1 -t , | 

5 

t. + T t.+T 


1+ 5 


In 


1 ' ' V+Tt-C') l+stv) 


n f 1 a° /V? &sisii 

71 


r. 


t. 


it l> lr i±T±lL 

- n Jf 
t f ) 

(7.2) 


dT : T-T' L 5 2 ( l: ).. e 2 ( l: .) “ 1 +E 2 ( T ') 


= S 1 + S 2 

where the first two terms are denoted by S-^ and the last one 
by S 2 - 

The shape of the path 5 (t) is shown in the diagram below $ 
and is derived in the text. 
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(*(t) ‘ % + 5 0 ) 2l“ 2 'VL < 5 < «. 


where ? + = 5 n ± "Ho 


(V.3) 


Putting (V. 3) in 


T 


:j o 

m|-2~ 

“(1 + 

m | pp” 

(5-^+5 

0 )| -(B) -- / 

0 


r. 

0 


T ° 

JdS In | 

2- n f 

0 

_2_ ( r 

2 C 0 U 

-ii+e 0 )| (B) +J a5 

in 

* 2 E 
^0 


’o 


(A) + f ° dC in |™(-l+p|-) | *(A) + Jd^n |^U 5 Q ) I U) 
o 2 ' 0 q 

(V-4) 


1+5 ‘ 


'fi ? Hr £ 

j o o 


5 + * ■ ’+ 

where A = -rr“ p in " 5 - = p - p 

51 -s’ 1+5 (iiA+5 rt ) --5' 


In 


i+(-n’+? 0 ) 


5- 


1+5* 


and B = — p In p* 

5 ^- r 1 + 5 * 


'o o 

BA “ 5 


1 + 5 


,2 


H- { p ! *» £ ) 

o ’° --- in 0 Z°L 


(t,;- 5 0 ) 2 - 6» 2 


1 + 5 


,2 
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She term proportional to In | [ in (7.4) is 

o ^o 5+ 

/ (A— B) d£ + / (A-B) + / (A-B) d£ . 

e- 0 2i ^ 


f+ s + i+c 2 

= / d? (-p- — y in ----- 

5. Sf-S 1+5 ‘ 


_£ 

5^-5 


c ‘ "2 1:: 


+ 5 


1 + C‘ 


(V. 5) is the same as exponent of the long, time term. 


ihe remaining terms in 3^ are. 


/ d£ lnj 


5 “P'+C 2 °o - 

— 1 • ( 3) - / d? • In |l-r^-~| • (3) 


5- 


o 


o 


*+ 


r -p'+l 

5 *o c 


5 --iA-e r 


f T _ I ■» ‘0 ’0 /»,\ , /> . _ ^ '0 •’O 

i d? In I™ *(n> -r J Cl? Ill — 


2n> 


5o 


I (A) 4 


2r) 0 e 

+ / d£ 111 I 1 - p— 1(A) + J d£ In I 
o ^o 2r)« 


5 “'i 0 “? 0 


(A) 


Collect midcle terms in both rows in (V. 6) 

9 -T 1 

2rj! - “ °o 

/ d£ In |1 - 27 - 1 (A) - c.5 in 

o ?0 0 


" 27 


( 3 ) 


¥e have, 


A-;-3 


A-B 


B 


A+B A-B 

2~ “ ~ F" 


So the above expression can be written as , 


/ 


° 0 d5 4,-"-ln |i-4 1 -i- / ° dC 7 =” In ( 

2 o 2 


25, 


2£ 0 + 5 


(7.5) 


(7.6) 


(7.7) 
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Bow take up the exprescion 


t - , . d £ ( u ) P 

j aT . aT . ^r T , 


5 ! 1 + 5 ^( t ' ) 


lil± 


In 


i r(r) 


Let us find the value of 


5( "0 



Divide % region into 4 parts 1, 2, 3 , 4 S 


(V.8) 


? (t) = 5° + and 5('C ! ) = 5 ° + t 1 

The expression (7.8) can he written as, 

( r ° -i- ci.%) a . - ( S ^ + a - r 1 ) a . 

i j J J r 


O 0 / \ 

= Si Cj - ( 3 =i + . 

First we check whether all the terms are indepenaent of t 


fhe only doubtful case is 1-3 9 3-4-° 

25, 


5° = ? 


55 


45, 


p CC— * — " r "* ) 

9 1 T 




? 
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If t is in region 1 and x 1 in region 3, 


t—t' = 


(t) 


5 s (~' ) “ 5e. 


o o 

^! a 3 5 3 a l 


5„ Sr 


4 5, 


( 4 - 6 ) = -2 


o 

5-jc 

’*'0 


o 


ax * dr ' 


85 


d5 ' ( t) • d5* ( x 1 ) . 


o 


Putting all tills in above ire see that r dependence cancels out. 


Similarly for 3-4 


How collect all the terms and put them in two classes I 
and II. Class I consists of those that have x,x' both in same 
region (1, S s 4) ana class II contains those that have r,r' in 
1-3, 1-4 or 3-4 . 


Terms in Class I 


(a) %,%' in region 1, 

2 5, 


£( x )- E + ~~~ X -4lLli _ 9, 


dr* = — d5 (t* ) 


25 

5(t)=- 5_ + 'C-'r 1 


o 


25 


- ( 5(t) - 5(t ! ) 


Then 


-1C = !° dc /V.^r . 

5 5 


1+ 5 2 


In 


5-5' 5 2 "?' 2 1+C 2 


(V.9) 


(b) Similarly, 

- txS ^ (3) = / as / 

0 o 


d5’ 


5 “5 


2 2 

5 -5' 


In 


It. 5 2 
l+5 ! 


(V.10) 
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(c) And -?eSp^ = / d£ I d 5 


5 + 5 + 


In (V.ll) 

. >~ c 4 l 


e; +5 _ e ++5 _ 


E ~ E’ p > p 

* * 5-5 1+5 


Terras in Class II; 


(a) TjjT* in regions 1-4 and 4-1 


25 


o 


5 (t) = 5 _ + pr— t, 
T o 

2 % 

g * ( t • ) = ^ + — -9- % , 


*. integrand is exactly sas 


,1(1) 

2 


ne as Sf'"' excepx region 


of 


. oil (1-4) , q II(4-l) _ 
^2 1 “2 


5 + 


5 ., 


. S^ J - VJ -^ r/ + = / 04 f dZ'(l) + f d5 / dC (I) 


s_ 5 + +« 


£ + +5 ~ 5 _ 


with. I = integrand In (V.9) . 


o !?+ 

= / d 5 J dg* 


5 __ £_+£__ - 5 * 2 


1 m i-±l! 


1 + 5 ’ 


2 * 


(b) t,t' in regions 1-5 and 5-1. 

2 5 Q 

5(t) = 5 + ~~£ t. 


“o 

45 


o 


5 ( T ' ) = 5 £_ + 'c . 

T o 


from page v 

A I(1 - 3) + si 1 ' 3 ” 1 ' = 4 / ««<•*) / ++C " «(*')<« + 

^ 851 5 0 “ T o T o 

4§ 

- s o 


85 o 5 


f_u\ P _ I 

H ; 2 5TT5 - 5( t ' ) +5 1 * w_\2 


ln l* 5 (t) 

5 (t)^- 5 (t ’) 2 1 + 5 2 ( t * ) 
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2 5+5 


+ 


J 


85? 6 


+ - o 45 25 45 

O • c / r , .. 


<i5(T) / d£ (t*) (^5_+~- 2 ^)t-' ? -* (*—•) 
5 o X “o 


5TtT^P5 


5+5 


1-n 2-±L±Ei. 

— — *- J - o 


X 


: ( + c(t) 2 -{(-e') 


= / cLS(-r) /' ac(t') -In ^ 

F o -5lty-5i.T j+5__ UF ^,, 


4 5 n 4 5 n 

{ 25 + ~T -35 - — > 

T o “ T o 

5.+5 

P . 


1+5^(t’ ) 


4- 


r / d5(T) / d C ( T * ) 


1 


i «r^TT0c. -^77^7) 


m ^cfh> , + / 0 ds ( T) /' : " E_ d ; ( T 1 ) 

1+ 5 1 ^ ( t 1 ) 5 o 


— — — X — 2 n XtiL XjeI 

5' 2 (t)-5 2 ('c } ) 1+5 2 (t') 

(Y.12) 


(c) t,t' in regions 3-4 and 4-3= 

5. 5.J+5 

sf I(3 “*> + sj^-V . / M f - dC ^ 


In 


5+ 5X5 


-c / + «7 


d5' 


V~£_ 0 


5, + 5 o 

T — 


2 £ — £ J ~4- £" 


l+_5_t 

1+5' 


-JL-tt m 

5-5^ 1+5 r 


(V.13) 


Anci that completes the list. Now we have to evaluate them. 

Most of the integrals are of the form 

b d ln(l+5 2 /H- 5' 2 ) 

/ d5 / d5' 
a c V 


? 2 

-4. J-l ^ 


let 1(a) = / d£ /d£' ixLi±XX/Xh a XI 


2 ,2 

5 "5 
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dl 

da 


= Jd£ fdV 



1+aS 2 


Jj_ 2 

1+a £ 



= / 

a 


1+a V 


* ! 


a 


_d 

_ 

1+a €' “ 


(tan 


■a 


|V o — 


-1 


—1 —1 
;aii ad-tan ac) 


Required integral = l(±) = ±(l) - l(u) 


da i 


-1 


/ WWU f t _L_ -T 

— ('can ccb-'c; 
^ n 
o 


"ca) ( tan"+" , ~' ! ' 


o:a~‘can ac, 


We next evaluate t q in to nos of moael parameters. the only 
t 0 -dependent terms in the short-range term is a In Jt q J . 
'ihis has to bo added to the contribution from the adiabatic 
plus Gaussian factor and then minimised. 
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Prom (page 69 ) A = 


i f dr rf • tan " 1 UiUKSnl 
^0 L 2 a 


+ 2TY r (T)TarTT"5) 111 {1 + (---^ ~ a “^~) 2 }] • Y ’ ( r) +crX( r) ) 

(V. 14 ) 


and the exponent of the G-auss ian factor is 


Y = -/' + 


2 ( _, -~2 r 

U u 


p /£lil + Ifiii) 


- -/aa^ - 


(V.15) 


Adaing these two contributions , 


r ,P 

Z A +y =- / ar V(X,Y') 
cr 0 0 


where, V(X,Y') = ^ 


2 <Y '^li + i 1 / 


ede 


U 


a — ii ( e — 1 ' — oX) 


(V.16) 

(V.17) 


She time arguments in X and Y' are understood; the last 
integral in (V.17) is the same as RHS of (V. 14) . The integral 
in (V.16) is, 


/ P dr V( X, Y' ) = rV(X, Y f ) | ^ - Jd-r-T dV 
o o 


Cl'U 


dV 

pv^ 0 ,x 0 ; -j ut-x- e 


= 6Y(X„,Y') -/ 


( V. 18) 


(X Q and Yq are the values of X, Y’ at either minimum of 
free energy.) Since g~- is non vanishing only if or is 

is nonvanishing, i.e., only during hops, (V.18) is 


/ 


dr V(X, Y’ ) - pV r 


h +T 0 

-2 n J 


dr • t* 


dV 

dr 


(V.19) 
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since there are 2n hops, the contribution from each hop being 
the same. Hie integral on the RHS of (V.lfc) can also be 
written as 


/ = -h- /U-ffcC 

O —A, 


(¥. 20 ) 


o 


substituting for time dependence of X from (¥.16). 
Again from ( ¥. 14 ) and ( ¥. 15 ) 

?(X,Y f ) = ~~ - " + ~ I I - , 

a -- jj ) + A 


X > £, + 


U 

1^2' 


_ (z ~ e a ) l t i |° sfe , a r° ete. fxl < e +- 

*-£ E 2 +1 2 * I (c-2I) 2 +4 2 ’ ' ' ^ 2 


u 


u 

2 ” 
(V.21) 


So 

that 

,u, .... 

,U v 

d¥ 

2X 

~ i [ tan" 1 .tan- 

71 u A 

-i fd + r*i i 

A “ 

dX 

~ U 


2£ d 
~ U 

. 2 r% j. ”1 2Xi | 

+ tc *-2 “ Lan a 3 s 1 

< £ d + 2 


- + 5 

•d 2 


(¥ 122 ) 


The integral on the r.h.s, of (¥.20) now becomes, since it 
is symmetric in X, 


x 


2 f 0 X p4 dX 
o UA 


d¥ 


,U T u 

E d + 2 2 

2 / X -~t dX + 2 f X ™ dX 

o ,U ^ 

£ d"‘~2 


(V. 23) 


Prom (¥.21) and (¥.22), the first integral in (¥.23) is. 


s 

/ v d¥ , y 

' A dX dX n 

o o 


/ £q 2 { ( l+-~) X - §X tan" 1 — } dX 


71’ 


U 2 n , 2e dx3 A ^ £ d + 2^ 1 , Un 2, a2 , -1 2£ d +U 

8 l” 1 ” ~U + x n { (£ d + 2 } + 4“ }tan ” A 
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which reduces, for a highly asymmetric raoael (U+ 2 £ d >>a ) 
U 


£ d + P 

t f. dV „ 
} a ^7 cu. 'v 


£ d U e 


3 


ax 


a 2 f_ £ d. AU _ A?, 

— 4 ‘ ti £ d 2'jx"" 4ic " 87 c 


(7.24) 


ihe second i nt e gr al in ( V . 2 3 ) 

U/2 d 7 > U/2 2 - >'2 e ^ + ~ r ^ 

X ~~ dX = j [—• - ~ {tan 
- , TT/o dX ,n/n d 7C 

e rl +U/2 e^+U /2 


/ 


■cT2 

'a 


+ 


'd 


•d 


U 


+ tan 1 JlJL}] dX 


.2 2 

3U 1 ^d 1 1 £ d“3’ ,c- d 


2 I ej{ ^4le„| n + |d 2 } + — li 


7C 


Afe d | ^ ( 2 I £ d I ) A(U-|e^|) 

It 


4|e d |(U _| E d| ) 


+ 


(7.25) 


for the highly asymmetric case. 

Adding (7.24) and (V. 25), we get for the case U>>|e d J»A 

,3 |e,|U 


/ X 


dV 

dl -UA= "5 


dX = 1 


~d 1 a_ 
U “4 


(V. 26) 


Putting (V.26) into (7.23) into (7.20) into (7.16), and putting 
the minimum condition, we get, 

4 


T o^ 


I I / £ * e d ' \ 

l £ d' ^ “ 3 tj 2 ' 


2 



APPEIjuIX VI 


D 


THE INTEGRAL I 


/. } - f 

v j ~ 1 x + a 


-jj 


Tills integral has be in. av... lusted below with tha help of 
contour integration. Hare f(x) is the Fornii function — 


ICO = f V tt— 

-L> (l+e 1 }(x+a) 

Wo me.-n to take the orincinal value of -™- 

x+. ■ 

integration variable in above- 


l+B*' 


ranging the 


I (a) = / 


PD 


dx 


-pD (l+e x ) (x+pa) -p jj (l+e^)(x+pa) 


= / 


dx 


pi 


dx 


(l+e x ) (x+pa) 



o 


dx 

Z i x 
1+e 


[i 


i 


x+pa 



o l+e x x+pa 


x u _, „ _ x-pa 


+-L-i 


/ 


-f)U 


O 


dx 

s-pa 



j ^ dx 

o (l+e X )(x+pa) 





(AVI.l) 


(.+VI.1I) 

( i'-.VI • 2) 


where the contour for integration for z runs fro m origiLn 
along real axis to pD. 

The contribution from above is the same as from the contour 
shown in the upper half plane since no singularity of the 
integrand falls within it. So (AVI. 2) can be written a s 

/ n — „ + f hx_ _ _ 4 x27tl 0 /e^£H 

o (l+e 1 ^) (y-ipa) o (l-e x ) (x+irc+pa) ^ pD 


(AVI. 3) 




24. Contours for the mTeoraftoh 
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Ihe third term comas from the integral around in and the 
fourth one is error due to extending the integral unto 00 , 

assuming | a | << lu 

In a similar way, the contribution from the Inst term of ( AVI . 1) 
can be shown, taking along the contour shown in the lower half 
plane, to be 

2rcixv a , 

+ 0 (2—) 


r w + r ta__ fit + o ( *± D - 

o (l+e iy ) (y+ipa) o (l-~e^) ( x-irc-pa) pD 


(AVI *4) 


adding (AVI- 3) and (AVI. 4) we get 


r % _.dy„ „ M + p t 

i y*~ipa irt+pa J 


xdx 


Tt+pj 


o (1-e )’ {x +(7t-ipa) *■> 


(AVI. 5) 


The last integral in (AVI. 5} can be expressed in terms of 
digamma functions 37 and we have the expression as, 


la ( 4|a±a) 


X'JX 


In 


( =k ~' r 4_ — i — « 4- ]h ( 1* 4- ) 

' O r O/rr-? * ^ Q. o 1 y \0 ^ O' iri- / 


i7i4-pa ^2 2n± ' 1 ^ + |a ‘ Y v 2 2%± J 

i % 


- ln <rff + ♦ + Ilf 


m (- 2 -f a ) - 2* (j>f> - *(|fr> 


(AVI. 6) 


Collecting (AVI .l*} and (AVI.6), we finally have, taking real 

parts since only they will survive, 

\ 

1(a) = ln M + ln r |- + 2He * (ff-) - * (fff 
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or 1(d) - In (f§) + 2Re ♦ (jfj) - Re * (^g.) (A. 71.7) 


where we have put 

k o 

— ; 



The low :.nd 

high 

t ar.ro 9 :^ turo 

nsions of (AVI 

.7) are, 

from small and la: 

rge 

y expressions of 4 (ig) 


1(a) .= In 

Jj-J. 

L> 

- | (|A 2 , 

<< a 


1(a) = In 

<*3> 

+ i- (5) (fj~) 

2 

? tcT >>a 


where q = ~™ 

2e r 

» r 

= .57721 56649, 

(3) = 1.08232 

32337. 
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